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Interpolation Properties of Orthogonal. Sets of Solutions 
of Differential Equations.” 


By O. D. KeEttoae. 





As the present paper is of the nature of a continuation of the one pub- 
lished in the last number of this Journal, an extended introduction may be 
dispensed with. Suffice it to say that we shall here be concerned with the 
problem of extending to the orthogonal function sets arising from ordinary 
differential equations of second order, the properties there derived for sets 
arising from integral equations. 


1. The Differential Equation and Boundary Conditions. 


We shall be concerned with differential equations of the form 


£ (p(2) 4) + (a(2) +ar(2))9=0, 


p>0O for 0<a<1, r>0 for 0<27<1, (1) 
with the general homogeneous self-adjoint boundary conditions 
ay (1) + by’ (1) =cy(0)+dy’(0), a’y(1) +0’y’(1) =e¢’y(0) +d’y’(0). (2) 


If p(x) vanishes at one of the end points of the interval (0,1), but in 
such a way that the differential equation has a solution which remains finite. 
and different from zero, with a finite derivative in the neighborhood of this 
point, one of the boundary conditions (2) is to be replaced by the condition: 
y(x) remains finite and different from 0 in the neighborhood of this point. If 
p(x) vanishes at both end points, both equations (2) are to be replaced by this 
demand for both end points. Examples of these cases are: (a) the functions 
J (ax) Vx, Jo(a,x) Vx, ...., Where a), a,,.... are the successive roots of the 
Besse] function of zero order J)(x); (b) the Legendre polynomials on the 
interval (—1, 1). 





* Presented to the American Mathematical Society, December 1, 1917. 
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The existence of solutions of the problem (1) and (2), and their oscilla- 
tion properties have been extensively studied,* and we shall make use of some 
of the more common ones. Our interest will center rather in the interpolation 
properties, and some of the consequences of the property (D),+ which it will 
be our object to establish. 

By “harmonics” we shall understand solutions of the differential equation 
and boundary conditions. The corresponding “frequencies” are the values of 
4 for which these solutions are possible. There is an infinite sequence of 
frequencies, without finite limit point, and bounded below. We shall think 
of them written in ascending order of magnitude, a frequency being written 
twice if it corresponds to two harmonics. The harmonics are orthogonal; ¢ 
even two corresponding to the same frequency may be made orthogonal by a 
proper choice of the integration constants entering them. No harmonic van- 
ishes simultaneously with its derivative. The zeros of any harmonic separate 
those of any harmonic with the same or lower frequency. In addition to these 
known facts, we assume explicitly (1) the continuity of the coefficients of the 
differential equation and of the derivatives involved, (2) that the boundary 
conditions are such that for any pair of harmonics F 


\ | (0), 9; (0) $:(1), (1) 

“iad Bet iia igi (3) 
| 94(0), ;(0) 9)(1), 95(1) 

and (3), that the i-th harmonic has exactly 1 zeros in the interior of (0,1) for 

all 7. A harmonic is “even” or “odd” according as the number of interior 





= p(1) 





zeros is even or odd. 
2. Types of Boundary Conditions. 


For what follows it will be convenient to separate the boundary conditions 
into types according to the behavior of the corresponding harmonics at the 
end points of (0,1). Consider first the case in which the determinant of 
coefficients in (2) ab’—a’b vanishes. There will then be a relation of the form 
a’’p;(0) +b’9;(0) =0. From this it follows that the function 


M (x) =p (x) (9; (x) o,(x) —@, (x); (2) ] 





* See Bocher, “ Lecons sur les méthods de Sturm,” Borel monographs, Paris, 1917. Rich indications 
as to the literature are found in the footnotes. 

f AMERICAN JOURNAL OF MATHEMATICS, Vol. XXXVIII, No. 1 (1916), particularly (1), (2), I, II, 
IV, V, VI and V’. Sturm (Liouville’s Journal, Vol. I, p. 433), and Liouville (same Journal, Vol. I, p. 269) 
studied the interpolation problem, and the method of the latter has suggested that used here. But the 
only existing results appear to be confined to boundary conditions of Class I below, which does not include 
the periodic and other interesting cases. 

¢ Orthogonal here in the sense that So (x) p;(a)r(w)da=0, ix}. 
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vanishes for x=0 for every 7 and j, and in consequence, by (3), also for «=1. 
We are thus led to the first class of boundary conditions, in which is also 
included the cases in which p(x) vanishes at one or both end points. 

Cuass I. M(0)=M(1)=0. One boundary condition involves only one 
end point, and the other, only the other. Consider first the case in which p (2) 
vanishes at an end point. Then under our assumption, all harmonics are 
different from 0 at that end point. Consider an end point at which p(z) does 
not vanish, say x=0. Then, as M(0)=0, $;(0)9;(0)—9,(0)9;(0) =0, and 
since no harmonic vanishes simultaneously with its derivative, it follows that 
if any harmonic vanishes for x=0, all do. We conclude for Class I: If any 
harmonic vanishes at an end point, all vanish there. 

Supposing the determinant cd’—c’d=0 leads to the same class. We may 
therefore assume that both this and ab’—a’b are different from zero, and that 
p(0)>0 and p(1)>0. The conditions (2) may now be given the form 


y (1) =ay(0) +6y'(9), (4,) 
y’ (1) =cy(0) +dy’(0), (42) 
or y(0)= dy(1)—by'(1), (5;) 
y' (0) =—cy(1) +ay’(1), (52) 


where, by (3) and (4), A=ad—bc=p(0)/p(1)>0. These conditions charac- 
terize the second class of cases. 

Cuass IT. M(0)=M(1)$0. The boundary conditions both involve both 
end points. M(0) and M(1) may both vanish, but not as a consequence of the 
boundary conditions. There will be three types to consider in this class. 

(a) b=0. Then ad>0, and (4,) shows that if 9;(%) vanishes at either 
end point, it does at the other. The sign of d now becomes important. 

(a,) b=0,a>0,d>0. If $,(x) vanishes at the end points, (4,) shows 
that $;(0) and ¢;(1) have the same signs, so that $;(7) must have an odd 
number of interior zeros. But (4,) shows the converse. Hence, for this case: 
All odd harmonics vanish, and all even harmonics are different from zero at 
both end points. An example of this case is y’+Ay=0, y(1)=y(0), y’(1) 
=y’(0), with the harmonics 1, sin 272, cos 27a, sin 47x, cos 4mm, ..... 

(a,) b=0,a<0,d<0. By similar reasoning, we conclude for this case: 
All even harmonics vanish, and all odd harmonics are different from zero at 
both end points. An example is y”+Ay=0, y(1)=—y(0), y’(1)=—y’(0), 
with the harmonics sin 2%, cos mx, sin 37x, cos 37%,..... 

(b) b>0. Since a harmonic will not vanish simultaneously with its 
derivative, (4,) shows that no harmonic vanishes at both end points. If 
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$;(0) =0, (4,) shows that $;(1) and $;(0) have the same signs, and the har- 
monic must be even. The same conclusion follows from (5,) when ¢,(1) =0. 
Hence in this case: No odd harmonic vanishes at either end point, and no har- 
monic vanishes at both. An example is y”’+A’y=0, my(1)=2y’(0), 2y’(1) 


=—zy(0). The harmonics are, for 1=2, $;(”) =sin a(o—1) += COS A,&, 
where A; is the i-th positive root of 4nA=(2’+4A’) sin A; while 2,=”,=7/2. 


For these frequencies we have the solutions A cos 3 za+B sin 3 x. It will be 


seen if A and B are so chosen as to give an odd harmonic, this harmonic will 
be different from zero at both end points. All later harmonics are different 
from zero at both ends. 

(c) b<0. By similar reasoning we conclude: No even harmonic vanishes 
at either end point, and no harmonic does at both. 


3. Lemmas. 

We shall need the following: 

1. If two harmonics have the same frequency, the one with the greater 
number of interior zeros is different from zero at both end points. For the 
zeros of the two separate each other. 

2. Let y(x%) =e.) (x) +0,9,(7) +....+0¢,,(%) be a linear combination 
of harmonics, in which at least one of the constants c; before the last is not 0. 
If, for c,=0, y(x) has p interior roots, then for some c,>0, tt will have at 
least p interior roots. If for c,=0 it has an infinite number of interior roots, 
then for some c,0, it will have at least N interior roots, N being any positive 
integer. The proof is easily supplied by considering the limit of y(x) as 
c,--0, the approach of y(z) to its limit being uniform. We shall reserve the 
notation y(x) for functions of the above form. 

3. Given a function y(x), ¢,—0, c, remaining fixed, we can choose 
Coy Cry ose +y Cay SO Small, not necessarily zero, that y(x) has the same number 
of interior zeros as $,(x) (1. €., exactly n), provided $,(x) does not vanish at 
an end point at which an earlier harmonic is different from zero. And in any 
case, y(z) has not more than one additional interior zero corresponding to a 
terminal zero of $,(x).. The proof depends upon the facts that $,(”) does not 
vanish simultaneously with its derivative, and that all the harmonics are finite, 
with finite derivatives. 





gta 2 
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4. Establishment of the Property (D). 


The results and treatment differ for the different types of boundary con- 
ditions enumerated. But the general argument is this. With a solution (2) 
of the differential equation, corresponding to boundary conditions to be 
assigned later, and not vanishing in the interior of (0,1), we form the differ- 
ential operators to be applied to y(z) : 


Ny (2) =p(2) [y'(@)V(@) —y(a)¥(#)], and Ly(2)=2Ny(z). (6) 


It will first be shown that repeated application of the operator Z reduces 
relatively the earlier coefficients of y(x) to any desired degree, and thus, by 
the third lemma, leads to a function with no more than m interior roots. It 
will then be shown that the application of the operator Z does not diminish 
the number of interior roots. The result, with the modifications noted, will be 
that no y(z) has more than » interior roots, and accordingly, that the deter- 
minant of the equations y(x)=0, y(a,)=0,...., y(%,)=0 does not vanish 
for any set of arguments x;, no two of which are equal. The property (D) 
will follow (see this Journal, Vol. XL, p. 146, footnote). 

Note first, however, that D,.(x,) >0, since the first harmonic has no 
interior zero. Second, that in case 2,=A,, D,(%, 2) >0. For the two har- 
monics, belonging to the same frequency, will have a non-vanishing Wronskian, 
from which fact the desired result may be inferred. We may therefore sup- 
pose in the following, n21, or, in case A,=A”), n22. 

The choice of ¥(2) in the operators (6) may vary with the function y(z) 
to which it is to be applied, and the corresponding parameter value 4* for 
which y(a) is a solution of the differential equation (1) will vary with it. We 
shall assume, however, that there are two constants, x and p, such that always 


Ata Agta 
. 1 3 -. | (7) 


The assumption will be considered below as use is made of it. Using the 
differential equation, it will be found that 


Ly (x) =p (a) d(x) r(x) [eo(A* —Ao) Go (%) +61 (A* —A1) oi (@) 
+... +c, (A*—A,)>,(2) ]. 
Dividing, on the assumption c, 0, by p(x)P(x) r(x) (A*—A,) we arrive at a 
new function 








x <A* <p where p< , or if A=A, p< 


uy (2) =65(—*) (2) + (% 


—) 1 (2) 


A*—2,, 
A*—A,_1 
eee bey a (M4) (2) +eaba(2). 
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* 
Under the assumption (7) the quantities = 





) are all less in absolute 


value than certain quantities independent of A* and less than 1, an exception 
arising only in the case of the last of them when 4,=A,_,. Hence, by repeat- 
ing the application of Z and the division, we may reduce all the coefficients of 
y(x) except the last (or, in case 2,=4A,_,, the last two) to be less than any 
assigned positive quantity. This completes our first step. 

The next step is to study how this process affects the number of interior 
zeros in the various cases enumerated. 

Cuass I. Choose ¥(z)=@(x). Then ~*=A,, and (7) is satisfied. If 
y(x) has p interior zeros, Ny(x) will be seen to vanish between each pair, and 
thus have at least py—1 interior zeros. But since in this case M(x) vanishes 
at both end points for all pairs of harmonics, and as y(z) is linear in the har- 
monics, it follows that Ny(z) also vanishes at both end points, and hence 
Ly(a) has at least p interior zeros. Thus no interior zeros are lost by the 
repeated application of Z. Lemma 2 permits us to assume c,>0, and then 
Lemma 3 shows that no y(z) has more than n interior zeros, since by § 2, no 
harmonics vanish at an end point under the present boundary conditions 
unless all do. And a case of double frequency does not arise.* Hence, in 
problems of Class I, the property (D) holds. 

Cuass IT, (a,). b=0,a>0, and y(1)=ay(0), for y(zx) is a linear homo- 
geneous combination of harmonics, and hence satisfies the same boundary 
conditions. Hence y(0) and y(1), if not both zero, have the same signs. We 
again use (x) =,(2), and assume that the last harmonic in y(x), $,(%) is 
even. If y(x) has an odd number, 2p—1, of interior zeros, it also has two at 
the end points, and Ny(z) has therefore at least 2p interior zeros, and Ly(z) 
at least 2p—1. Thus, if the application of Z to y(x) changes this function 
from one with an odd number of interior zeros to a new one with an odd 
number, none are lost. If it changes y(x) to a function with an even number 
of interior zeros, at least one is gained. If y(x) has an even number, 2p, of 
interior zeros, Ny(x) has at least 2p—1, but (3) shows that it then either has 
an additional interior zero, or else two at the end points, and in either case 
Ly(ax) has at least 2p—1 interior zeros. If Ly(x) has an even number of 
interior zeros, it must be at least 2p, and if odd, at least 2p—1. But ulti- 
mately, by Lemma 3, the repeated application of L leads to a function with an 
even number of interior zeros, since we have supposed this to be the case with 





* For two harmonics of the same frequency, M(#)=const. This constant is not zero if the har- 
monies are independent. But it must be 0 in the problems of Class I. 
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,(”). Hence no zeros have been lost, and we conclude: The harmonics of 
Case II (a,) have the property (D) for every even n. 

The conclusion is still justified in case ~,_,=A,. Here we must study 
C,-1P,-1(%) +¢,9, (2). If this has an even number of interior zeros, their 
number can not exceed n, since they are separated by the zeros of $,_,(”), and 
C,-1Pn_1(%) +0¢,,(2) will, by the boundary conditions, be different from 0 at 
the end points. The conclusion thus subsists in this case. If c,_1$,1(x) 
+c,,(x) has an odd number of interior zeros, the repeated application of LZ 
to y(x) may have reduced the number of interior zeros of y(x) by one, but by 
no more, and ¢,_,9,-1(%) +¢,,(%) has n—1 interior zeros, and none at the end 
points, since by § 2, @,_,(”) vanishes at both end points, and ¢,(x) at neither, 
and c,+—0. The application of the lemma 3 then shows that y(x) can not 
have more than v interior zeros, and the conclusion is established also in the 
case A,_,=A,. 

Cuass IT, (a,). b=0,a<0, and y(1)=ay(0), so that if y(1) and y(0) 
are different from zero, they have opposite signs. We now operate with Z in 
which (a) is chosen subject to the condition that y’(x)(x2)—y(a)V' (a2) 
vanishes at both end points. By § 2, @)(x”), being an even harmonic, vanishes 
at both end points, and as (x) is different from zero on the interior of (0,1), 
it follows that ~* <”~,. Hence the assumption (7) is justified as far as the 
upper limit on 4* is concerned. The existence of the lower limit will be 
generally established at the end of the paper.. We start with a function y(z) 
ending in an odd harmonic. The operator N does not reduce the number of 
interior zeros by more than one, while with the present choice of Y(x), Ny() 
has zeros at the end points. So the operator Z does not diminish the number 
of interior roots. As the odd harmonic ¢,(x), by (2, is different from 0 at 
the end points, Lemma 38 gives the result: The harmonics of Case II (a,) have 
the property (D) for every odd n. 

Cases II, (b) and (c). Here, again, we choose (x) so as to make Ny(x) 
vanish at the end points. It is not, in these cases, apparent that the upper 
limit for ~* in (7) holds. However, since ¥() does not vanish within (0,1), 
and as @,(”) vanishes twice, it is clear that A4*<A,. As the frequencies have 
no finite limit point, there must be a positive integer n, and a number , inde- 


pendent of 4*, such that A*<p < om, The least value of » for which these 


~_ 


inequalities are possible,* we will call n’. The type of argument used in the 





* This least value of can often be easily determined in a particular problem, it being helpful to 
observe that the value of \ for which the differential equation has a solution vanishing at the end points 
of the interval (0, 1) but not in the interior, is an upper bound for A*. In the example in §9 under Case 
II (b), (7) is fulfilled without alteration. 
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previous cases then permits us to conclude: The harmonics of Case II (b) 
have the property (D) for every odd n2n’, and, provided none of the har- 
monics after dy_,(”) vanishes at either end point, for every n2n’, -Also, the 
harmonics of Case II (c) have the property (D) for every even n=n’, and, 
provided none of the harmonics after $,-;(x) vanishes at either end point, 
for every n=n’. iil 

5. Existence of a Lower Bound for 2*. 


In the problems of Class II, from (a,) inclusive, on, it was assumed that 
a lower bound for the values of 4* existed. We proceed to show that such is 
the case. It will be recalled that ~(a2) was to be chosen so that the ratios 
Wy’ (1)/Y(1) and ¥’(0)/(0) coincided with the corresponding ratios for y(z). 
Let us call these ratios 7 and & respectively. They are not independent, but 
are connected by the equation, obtained by dividing (4,) by (4;): 


n=(ce+dé)/(a+bé&), where ad—bc>0. (8) 


We shall now prove the theorem: 


There exists a number A, independent of & and x, such that the solution of 
the differential equation (1) which satisfies the boundary conditions 


w'(0)=§u(0),  w’(1)=nu(1), (9) 


subject to the relation (8), and which does not vanish in the interior of (0,1), 
belongs to a parameter value A* >A. 
We assume that (0) >0 and p(1)>0. The theorem is not needed in the 


problems of Class I. ; 

We first prove the lemma: For fixed & and decreasing yn, and for fixed yx 
and increasing &, the value of the parameter A corresponding to (9) increases. 
Let u,(z) and u(x) be two solutions of the differential equation satisfying the 
first equation (9). Denoting the corresponding parameter values by 2, and A,, 


we find the following identity: 
1 
p (1) (uy(1) w; (1) 1 (1) us(1)] = (Aa—Aa) fue (2) (2) 7 (@) de. 


We are, in this paragraph, concerned only with solutions of (1) which do not 
vanish in the interior of (0, 1), so that it is legitimate to consider u,(x) and 
U,(x) positive for 0<2<1, so that the integral is positive. Dividing by 
Uz(1)u,(1), we have 


p (1) bn (da) —n(%4)]= AE fun(a)(z)r(w)dz, (10) 














eS pa ae 
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which shows that decreasing y corresponds to increasing A. The second part 
of the lemma is proved similarly. 

We next find the relation between & and 7 when these are the values of 
the ratios wu’(0)/u(0) and w’(1)/u(1) corresponding to any fixed a. Todo 
this we take a fundamental solution set for z=0, u(x) and u,(x), for which 
u,(0)=1, 4,(0)=0, u.(0)=0, w,(0)=1, in terms of which any solution u(x) 
for the same value of 4 may be expressed in the form u(x) =¢, u(x”) +6, U2 (2). 
Then £=c,/ce,, and y=[c,u,(1) +¢.u,(1) ]/ [eu (1) +¢,u,(1)]. Eliminating 
the constants, we find the relation 

_ uy (1) +§us(1) 
”= an (1) + E(1)’ sie 


for which we have the easily derived relation 
Uy (1) ue(1) —u (1) us (1) =p (0)/p (1) >0. 


This shows that the hyperbola which is the Cartesian representation of the 
equation (11) is an always rising one. For various values of A, this hyperbola 
varies in position and size, but if for a finite value A of A, its center lies above 
and to the left of the hyperbola (8), also an always rising hyperbola, the 
lemma shows that the solution of the differential equation and the boundary 
conditions (9) subject to (8) will always correspond to a value of A>A, and 
the theorem will be proved. 





That for some finite A the center of the hyperbola (11) lies above and to 
the left of the upper left-hand branch of the hyperbola (8), will be evident 
from the fact which we shall next prove: given any finite positive quantities 
h and k, a finite A always exists for which the coordinates of the center of the 
hyperbola (11) satisfy the inequalities 





Uy (1)/u,(1) >h, ug(1)/u,(1) >k. (12) 
; . U(l)  ¢’  p(0) ‘att ite ; ; 
The identity u, (1) a paie dx, the derivation of which gives no 


trouble if the initial values of the fundamental solution set are kept in mind, 
shows that the ratio u,(1)/u,(1) can be made as great as we wish if u,(x) can 
be made as great as we please at all interior points. But this fact can be 
inferred from the identity 


d (402 2} -_ A,—A 
dar U; (2%, A,) p(x) uz (x, ay) 


This shows that u,(”,”)/u,(x,”%,) imereases with w if A<A,. This ratio 
30 





Xx f(a, A) Uy (2, Ay) r(x) da. 
0 
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approaches 1 as x=0. Hence always, u,(%,A) >u,(%,a,) for 0<a<1, and 
using the law of the mean, o we a j r(v)du>(A,—A) X a posi- 
tive function of x independent of ~. Hence the ratio, and therefore u,(2, ”) 
can be made as large as we please, and the first inequality (12) is established. 

To prove the second, we start by noticing that for sufficiently small 4, 
u;(z)>0 for 0<x2<1. This is seen from the equation p(x)u;(x)=p(0) 





+f (—ar—q) ude, in which the integral is positive if A<—max. g/min. r. 
0 
Multiplying the differential equation for u, by 2p(x”)u;(x) and integrating, we 


1 1 ee ee 
have [pus(2)}*| =—2f p(Ar+q) ugugdv=—p(ar+q)uz(1), the bar denoting 
0 0 


a mean value. Hence ui(1)/uy(1) > = V(—A) min. pr—max. pq, which 


can be made as large as we please by sufticiently diminishing 2. Thus the 
second inequality (12) is proved. 





CotumBIA, Mo., October 20, 1917. 




















Directed Integration. 


By H. B. Pumuies. 





In case of an integral along a curve, 
fPdx+Qdy=Lim XPAz+ QAy, 


the increments Az and Ay may be positive or negative according as x and y 
are increasing or decreasing. In case of a double integral, however, 


SSf (x, y)dady=Lim UYf (2, y) Ardy, 


the element AvAy is usually considered positive, or at least, invariable in sign. 
This introduces difficulties similar. to those which occur when we attempt to 
banish the minus sign from algebra or analytic geometry. For instance, in a 
change of variable, it is necessary to assume that the Jacobian has an 
invariable sign. 

Physicists avoid these difficulties by introducing a cosine which is positive 
or negative as required. Mathematicians accomplish the same result by 
making a change of variable, thus obtaining an element of integration which 
need not change sign. I wish to show in this paper how the algebraic sign 
can be directly attached to the element of integration, multiple integrals being 
treated in this respect like curvilinear integrals. The equations for change of 
variable and those connecting line, surface, and volume integrals, present 
themselves much more naturally in this form. In this discussion I shall not 
enter into questions of existence and convergence. These matters are treated 
in practically the same way whether the integral is directed or not. I shall 
also consider only two and three dimensions, although the extension to higher 
spaces is immediate. 

Directed Regions.—A surface is called one-sided if it is possible to pass 
from a point on one side of the surface to a point on the other without passing 
through the surface or across its border. If this is not possible, the surface 
is called two-sided. A simple one-sided surface can be formed by twisting a 
strip of paper through 180° and bringing its ends together. 











ar 
ry erecta 
x 12 


ay 
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If a surface is two-sided, one side can be considered positive, the other 
negative. We shall assign a direction or sense to a region on a two-sided 
surface by assigning a direction around its border. In a right-hand system 
the positive direction is usually chosen such that an observer on the positive 
side of the surface finds the region on his left when he moves in the positive 
direction along the border. It should be noted that this direction does not 
belong to the border, but to the region of which it is the border. Thus a given 
direction around a great circle of a sphere is positive for the hemisphere on 
one side, and negative for that on the other. 

Surface Integrals—Let x and y be one-valued and continuous functions 
defined at each point of a region I on a plane or two-sided surface. Divide I 
into elementary regions, or cells, by two sets of curves 


x=constant, y=constant. 


Any one of these cells whose boundary is a simple quadrilateral with two 
pairs of opposite sides belonging to the curves z, x+ Az and y, y+Ay will be 
called regular. Irregular cells may be bounded by less than four curves, or 
by four curves that are not of this simple type. 

In the definition of the integral only the regular cells will be used. Hence 
it is assumed that the irregular cells can be enclosed in a region or set of 
regions whose total area approaches zero when Az and Ay approach zero. 
This is certainly true of the curves ordinarily used in integration. It would 
not be true if the two systems of curves x=const. and y=const. were the same. 

Choose a direction around one of the quadrilaterals. Then for that 
quadrilateral we define AvAy as the product obtained by multiplying the 
increments of x and y which are found by passing around the quadrilateral in 
the chosen direction so as first to traverse a curve y=const., and then a curve 
x=const. The sign of AvAy is fixed for a given quadrilateral and a given 
direction around it. Thus in the quadrilateral ABCD, if AB and DC are por- 
tions of curves y=const., and AD and BC portions of curves x=const., we 
may take Az from A to B and Ay from B to C, or we may take Az from C to 
D and Ay from D to A. In the second case the signs of Az and Ay are both 
changed, and so AvAy is not changed. Similarly, the product AyAvw is 
obtained by traversing the quadrilateral in the same direction, first traversing 
a curve x=const. and then a curve y=const. Inspection of a figure will make 
it clear that one of the increments in AyAwz differs in sign from the corre- 
sponding increment in AvAy. Hence 


AyAr=— AxvdAy. (1) 
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Let the same direction be taken around all the quadrilaterals into which 
[is divided. Then, if (2, y) is any point in the quadrilateral to which AvAy 
belongs, we define the integral of f(z, y), 


SSf (a, y) dady, 


in the chosen direction over [ as the limit (if it exists) approached by the sum 
LU (a, y) AxAy 
when Az and Ay approach zero, the summation being for all the regular cells 
within. Similarly, 
SSf (a, y) dyda=Lim XX (2, y) AyAz. 
Therefore, by equation (1), 
SST (%, y)dydx=—JS ff (x, y)dxdy. (2) 


It should be noted that the order in which the differentials are written does 
not indicate an order of integration. In fact, no order of integration is con- 
sidered. The integrals are multiple, not iterated. . 

We have assumed that Az is determined along the curves y=const., and 
Ay along the curves x=const. It is a very important fact that one of these 
increments could be determined along a third set of curves w=const. Thus, 
if we resolve [ into cells by the curves x=const. and w=const., and in each 
quadrilateral determine Ay on x=const. as before, but Av on w=const., the 
value of the integral will not be changed, provided the total area of the 
irregular cells formed by the new curves has a zero limit. For TI can be 
resolved into strips between consecutive curves 7,x+Az. All the quadrilat- 
erals in a strip have the same Az. Also Ay is taken in both cases along the 
curves x=const. Hence, in the change assumed, the part of the summation 
belonging to this strip is affected only through the change in the distribution 
of the intervals Ay. This does not affect the limit. | 

Volume Integrals.—Triple integrals are defined in a similar way. Let 
x, y, 2 be one-valued and continuous functions defined at each point of a region 
I of space. Divide I into cells by means of three sets of surfaces x=const., 
y=const., and z=const. We shall call the cells regular which are bounded on 
opposite sides by three pairs of surfaces ~ and x+Az, y and y+Ay, 2 and 
ze+Az. We assume that when Az, Ay, and Az approach zero, the total volume 
of the irregular cells approaches zero. 

Let AB, BC, CD be consecutive edges of a cell, y and z being constant on 
AB, zand x on BC, x and yon CD. Let the outer surface of the cell be con- 
sidered positive. The path BCD determines a direction (positive or negative) 
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about the face of the cell in which B,C, D lie. I call this the direction, or 
sense, of the cell ABCD. 

For a given cell taken with an assigned direction, we define AvAyAz as 
the product obtained by multiplying the increments of x, y, 2 which are found 
by passing along consecutive edges of the cell in the assigned direction, first 
traversing an edge on which x alone varies, then one on which y alone 
varies, and finally one on which 2 alone varies. It is easy to verify that this 
fixes AvAyAz not only in magnitude, but also in sign when the cell is given 
and a direction assigned to it. Thus let the cell be a rectangular box 
ABCDEFG formed by the parallel rectangles ABCH and FEDG (x varying 
on AB, y on BC, and zon CD). We may take Az on AB, Ay on BC, and Az 
on CB, or we may take Az on EF, Ay on FG, and Az on GH. In the second 
case two signs (those of Aw and Az) are changed, and so AvAyAz is not 
changed. 

Similarly, the product AyAzAz is obtained by traversing the edges of the 
cell in the same direction, first traversing an edge on which y alone varies, 
then one on which 2 alone varies, and finally, one on which & alone varies. 
The other products of Av, Ay, Az are defined in a similar way. It is easy to 
verify that each inversion of the order of Az, Ay, Az introduces a negative 
sign in the result. Thus 

AxAyhz= — AxrAzAy= AzArdAy. (3) 

Let all the cells in a region [ be taken in the same direction. If (2, y, 2) 
is any point within the cell to which AvAyAz refers, the integral in the 
assigned direction over I, 

SSSf (x, y, 2) dadyde, 
is defined as the limit (if it exists) approached by 
LUX (x, y, 2) AvAyAz 
when Az, Ay, and Az approach zero, the summation being for all the regular 
cells within I’. Similarly, 
SSSf (x, y, 2) dydedx=Lim LUX (a, y, 2) AyAzAz, 


etc. Hence, from (3), 
SSIF (2, y, 2) dadyde=—fSfff(x, y, 2)dxdedy=ffff(x,y,2)dedady. (4) 
As in case of double integration, one of the sets of surfaces x=const., 
y=const., z=const., can be replaced by a third set of surfaces, provided the 
irregular cells thus introduced have a total volume that approaches zero in 
the limit. 
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Expression of an Integral in Terms of Boundary Values.—A simple 
integral is expressed in terms of its limits by the formula 


fdu=u 


Similar formulas apply to multiple integrals. 
Let u be a function of x and y. In the integral, 


Sfdudz, 


let du be taken along x=const. We may take dx along w=const. or y=const. 
as we choose. The integral of du along a curve x=const. from one intersec- 
tion with the border to another is given by equation (5). Hence, if we 
evaluate the integral by summing first with respect to u and then with respect 
to 2, we get 


“=u(b)—u(a). (5) 





Sfdudx=fuda. (6) 


The double integral is taken over a region I, the simple integral over the 
boundary of [. Since du and dz occur on consecutive sides of a quadrilateral 
in the double integral, the direction of integration around the boundary must 
be such that if BC belongs to the boundary and ABCD is a quadrilateral 
directed as in the double integral, then the integral along the boundary is in 
the direction BC. | 

Similarly, if u is a function of a, y, 2, 

Sffdudady = ffudady. (7) 

The triple integral is taken over a region [, the double integral over its 
boundary. The directions of integration are so related that if ABCD is a cell 
of the triple integral with face BCD in the boundary, then BCD gives in that 
face the direction of the double integral. 

Illustrations of these formulas are furnished by the theorems of Green, 
Stokes, and Gauss. Suppose, for example, P, Q, R are functions of 2, y, 2 on 


a two-sided surface. Then 
fPdx=ffdPda, 


the two integrals being taken over a region of the surface and around its 
boundary, respectively. Since dP is determined on the curves x=const., 


oP oP 
dP = = dy+ = dz, 


oP oP 
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In these integrals dy and dz are taken along x=const. In the combination 
dydz we may take dx along y=const. and in dedx along z=const. Similarly, 


sQay= ff 2 dandy + SP dedy, 


R 
NOM pcos Da a dz + dydz. 
Adding these equations and changing the sign each time we invert the order 
of differentials, we get 


P P R 
fPdz+Qdy+Rdz= (2 Oa me. oy tt dd y+(F— ay dz em —  \dede, 


which is Stokes’ Theorem. Since 


_ @P oP oP 


dP = ae a =| 


the equation 
siete 


suggests that symbolically 
f - = ~ dxde= 0. (8) 


This is in line with our definition since one set of curves does not give rise to 
any regular quadrilaterals, and so the summation from which the integral 
might be defined is zero. 

Change of Variable.—Let x, y be functions of u,v. Then 


Ox Ou 1, 
da= = du ut ard 


along any curve. Hence 


SSidedy = 752 dudy + fpf avdy. 


In these integrals dy is determined along the curves x=const. Since this is not 
significant we may determine dy in the first integral on the curves v=const., 
and in the second integral on the curves u=const. In the first case 


dy= £Y ae, 
and in the second 
dy= oY a, 
whence 
Lyme) Ox 0 Ox O Ox O 
SSidedy = p6t 5 Y dudo-+ 5 p¢3" SY dvdu= spr] 2 SY — 9% 4 audv. (9) 
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This result could be obtained by using the values 


Ox Ox 
dz= a, aut om dv, 
9 7, Py 


in fffdxdy, expanding, and assuming that the integrals containing dudu and 
dvdv are zero. 
In a similar way we show 
Ou dy de 
Ou du du 
Sf (fandydene ppt t ee S lain (10) 
dv dv dw 
Ox Oy de 


dw dw dw 








Equations (9) and (10) are valid whether the Jacobian has an invariable 
sign or not. Irregular cells will, however, usually occur in the neighborhood 
of a point where the Jacobian is zero. It may therefore be necessary that it 
be possible to enclose such points in an area in case of a double integral and 
in a volume in case of a triple integral which has a limit zero. 

While an integral is represented by a number, the summation process 
upon which the integral is based is as much geometric as arithmetic. Sym- 
bolically the integral is a function in which daxdy is equivalent to —dydx and 
dz’? to zero. This can be expressed by means of vectors. Yet dw, dy, dz are 
not vectors, and integration belongs no more to vector analysis than algebra 
does. They both belong to the larger field of quantities having sign, but not 
direction. , 
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P-way Determinants, with an Application to Transvectants. 


By Lerrne Hau Rice. 





In this paper an extended definition of a determinant is given which 
applies to determinants of more than three dimensions, and enables us to 
remove the restriction in Cayley’s law of multiplication and to set up a new 
case in Scott’s law of multiplication. New formulas are obtained for the 
known process of decomposition of a determinant into determinants of fewer 
dimensions, and a new process called crossed decomposition is described. 
Fresh light is thrown upon the function known as a “ determinant-permanent,” 
a limitation hitherto thought necessary being done away. Finally a generali- 
zation to p dimensions is made of Metzler’s theorem in two dimensions con- 
cerning a determinant, each of whose elements is the product of k factors. 

We lead up to these matters by a brief statement of the elementary 
theory of 3-way or cubic determinants and permanents. 


I. TuHuree-way DETERMINANTS AND PERMANENTS. 


3 


1. Definitions and Fundamental Properties.—Elements, n*? in number, 


can be arranged in a 3-way matrix of order n having n? rows, n” columns, and 
n* normals; all are called files. The matrix divides up into n strata, or layers 
that contain rows and columns and are pierced by normals; and into  row- 
normal layers pierced by columns; and into n column-normal layers pierced 
by rows; ail are called layers. In triple-index notation, a,, denotes the 
element in the 7-th stratum, the j-th row-normal layer, and the k-th column- 
normal layer. To represent 3-way matrices of successive orders, we write: 


Ay11 U2 Ayyg | Aer Mere Ae1g | A311 A312 Ag13 


A111 Ay12 | G11 G12 
’ Ay21 Ay22 Arg | A221 Aoe2 Moog | A321 Aza Agog |], ete. 


Ay21 Uy22 | %21 Aoo9 

















A131 A132 A133 | G31 Cege Aeg3 | A331 Ass2 Asse 

Two or more elements are conjunctive if no two of them lie in the same 
layer of any direction; conjunctive elements are perjunctive and form a 
transversal. When we speak of a determinant of a matrix and use the word 
“transversal,” we mean the product of these elements. The locant of an 
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element is the set of indices which locate it in the matrix. If the locants of a 
perjunctive set of elements be written. in a column, the three subcolumns, of » 
indices each, are called ranges, and the whole is called the locant of the set 
(transversal) ; @. g., Gio. A313 Aog,, ranges 132, 213, 231. The sign of a range is 
+ or — according as there is an even or an odd number of inversions of 
order in the range. 

The determinant of a 3-way matrix is the algebraic sum of its transver- 
sals, each having the sign which is the product of the signs of its second and 
third ranges. 

The permanent of a 3-way matrix is the sum of its transversals. 

The determinant and permanent being homogeneous linear functions of 
the elements of any layer, we have obvious theorems as to factors of layers, 
and as to separation into a sum of 3-way determinants or permanents when 
there are polynomial elements, and, conversely, as to addition of determinants 
and of permanents. An interchange of strata in a determinant, or of any 
parallel layers in a permanent, does not change its value; but an interchange 
of any two parallel layers other than strata, in a determinant, changes its 
sign. Hence, if two such layers are alike the determinant vanishes. Hence a 
multiple of such a layer may be added to any parallel layer without changing 
the value of the determinant. 

A minor is formed by striking out an equal number of layers of each 
direction. It is obvious what we mean by conjunctive minors and by per- 


_ junctive minors. <A perjunctive set of minors, formed into a product, with the 


proper sign, is equal to the sum of a certain number of terms of the determi- 
nant; the sign is the sign of that term of the determinant whose elements are 
the elements in the main diagonals of the minors. In the simple case of an 
element a,,, and its complementary minor, the sign is (—1)’**. In consequence, 
the Laplacean expansion of a determinant is formed by partitioning the matrix 
into two or more sets of parallel layers and forming all possible perjunctive 
sets of minors occupying the sets of layers, one minor in each set. 

2. Decomposition.—(i) A 3-way determinant A of order n can be decom- 
posed into the sum of »! 2-way determinants whose rows are rows of A. For, 
arranging n perjunctive rows of A in the order of the row-normal layers in 
which they lie, we see that the 2-way determinant 


Gey. Ayre +++ + Ain 
Aja Ojn29 eseee Ojr'2n 
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whose matrix they form consists of n! terms of A, and that the totality of 
such 2-way determinants consists of all the terms of A. They are called com- 
ponents of A. 

(ii) If the rows are arranged in the order of the strata in which they lie, 
then each determinant must have prefixed the sign of the j-range (denoted by 
+,) in the locant of the set of rows. | 

(iii) (iv) There will be two corresponding forms of decomposition into 
2-way determinants whose columns are columns of A. 

(v) We can also decompose A into an algebraic sum of 2-way permanents. 
Arrange » perjunctive normals in any order, and to the permanent whose 
matrix they form, prefix the sign of their locant; for, clearly, all of those »! 
terms of A which lie in a perjunctive set of normals have the same sign. 








Ayj1 Ayj2 eeese ay Qe eeee 
(11) £5] Goya Gaye. 43 (UT) | Qyer Ayre...» [5 
+ +% 
Ay Qo, eeee yj, Ae K 
(1V) Hg] Gop Qogy...|3 (V) ApH Dy jr pes Aggy ove |e 


A 3-way permanent may of course be decomposed by rows, columns, or 
normals, into a sum of 2-way permanents. 

3. Element-Multiplication.—(i) (Scott’st law of multiplication). The 
product of two 2-way determinants (or permanents) A and B of order n is 
expressible as a 3-way determinant (or permanent) C of order » wherein 


Cijp= Aj i - 


+ + 
Ay; %2 A104, Ay,D4 | Ae.) Ax Doo 


EXAMPLE: 


+ 
: by by; 
Doi doe 


From the prescription ¢,,=a,,b,, whatever be the value of », we see (i) 
noting the index i, that columns of A and columns of B are found in normals 
of C, (ii) noting the index j, that rows of A are found in columns of C and 
that elements of B run as factors through columns of C; and (ili) noting the 
index k, that rows of B are found in rows of C, and that elements of A run as 
factors through rows of C. 




















Az1Qo9 Ay2041 A420 2 | Age) A2qDo9 





* «| | ” means a determinant or permanent; «] 1” means a permanent. 

7 R. F. Scott, “On Cubic Determinants,” etc., Proc. London Math. Soc., Vol. XI (1879), p. 17, at p, 
23, paragraph 7. In paragraphs 8 and 9, Scott extends the rule so as to give the product of two deter- 
minants of p and g dimensions, respectively, in the form of a determinant of p-+-q—1 dimensions. In 
§7 of the present paper this rule is extended to determinants as defined in § 5. 
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This method of examining a prescription will often give at once a good 
idea of a matrix and may suggest the best way of dealing with it. 

In the present case we use decomposition (ii) of the preceding section 
and find that each component equals B multiplied by a term of A, giving 
AB=C. 

(ii) The product of a 2-way determinant A and a 2-way permanent P of 
order n is expressible as a 3-way determinant C wherein 


C je jx Dir s 
that is, a normal of C is a column of P with an a-factor; a column of C is a 
column of A with a p-factor; and a row of C is a row of A and a row of P. 


Using decomposition (v), we find that each component equals P multiplied by 
a term of A, and AP=C. 


4, File-Multiplication.—(Cayley’s* law of multiplication.) The product 
of a 3-way determinant A and a 2-way determinant B of order n is expressible 
as a 3-way determinant C of order n, wherein 


& 
= 
ll 

T Me 


: Oi D p13 


that is, the matrix of each stratum of C is precisely what would result from 
using the familiar process of multiplying together two 2-way determinants 
into a 2-way determinant, row into row; one of the determinants being always 
B, and the other being the determinant of the matrix of the stratum of A 
corresponding to that of C. In brief, we multiply B into the strata of A to 
form C. 

Using decomposition (i), we find that each component C,=—A,B, whence 
C= 48. ; 

II. DETERMINANTS AND PERMANENTS OF p DIMENSIONS. 


5. Definitions and Fundamental Properties —A p-way matrix (or a 
matrix of class p) of order » is formed of n? elements: 
|), 


I] Gajcg...-ctpl 
The matrix can be separated into n layers, in any one of the p directions; a 
layer is a (p—l)-way matrix of n?~’ elements. Common to two layers of 





* A. Cayley, “On the Theory of Determinants,” Trans. Cambridge Phil. Soc., Vol. VIII (1843), p. 
75; Coll. Math. Papers, Vol. I, p. 63. See §2 of that paper. Cayley extends the rule so as to give the 
product of two determinants of p and q dimensions, respectively, in the form of a determinant of p-+-q—2 
dimensions, but states that the rule is inapplicable when both p and q are odd; a restriction ignored by 
Scott in paragraph 9 of the paper cited in the first note, and by many others (see “Abrégé de la théorie 
des déterminants 4 n dimensions,” par Maurice Lecat; Gand, Hoste, 1911, pp. 55 et seg.). In §8 of the 
present paper this rule is extended to determinants as defined in §5, such determinants serving to remove 
the restriction. 
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different directions is a (p—2)-way sublayer; common to three, a (p—3)-way 
sublayer; and so on, until we have, common to p—1 layers of different direc- 
tions, a file of n elements, piercing the n parallel layers of the remaining 
direction. Parallel to this file are n?~'—1 other files, together with it con- 
taining all the elements of the matrix. Files of the last direction are rows. 

We speak of conjunctive and perjunctive elements, files, etc., of transver- 
sals, and of locants and ranges, as we did in the special case of a 3-way 
matrix. 

Hitherto a p-way determinant has been defined as the algebraic sum of 
the transversals of a p-way matrix, the sign of a transversal being determined 
by arranging its elements in such an order that the values of a fixed index 
shall read 1, 2,...., ”, and then taking the product of the signs of all the 
other ranges. It has then been shown that for a matrix of even class the same 
. determinant will result, whatever be the fixed index, but that for a matrix of 
odd class a different determinant, in general, will result from a different choice 
of the fixed index. It has also been shown that interchange of layers denoted 
by the fixed index in a determinant of odd class leaves the value of the deter- 
minant unchanged, while the interchange of two layers of any other direction 
in such a determinant, or the interchange of two layers of any direction in a 
determinant of even class, changes the sign of the determinant. 

It follows that without making the supposition that the elements of a 
transversal are first arranged in any particular order, we may say that the 
sign is the product (i) of the signs of all the ranges, in a determinant of even 
class, and (ii) of the signs of all but one (a fixed one) of the ranges, in a 
determinant of odd class. 

From this point of view, we now generalize the definition of a determinant. 
We shall call an index (or range, or direction, or file, or set of layers) signant 
or nonsignant according as we do or do not take the order therein into account 
in fixing the sign of aterm. Ina 2-way determinant both indices are signant ; 
in a 2-way permanent, both indices are nonsignant. In a 3-way determinant, 
two indices are signant. Passing to matrices of more dimensions than three, 
we see that it is possible not only to have signant all the indices if the 
class is even, and all but one if the class is odd, but to have signant a less 
number in either case, provided only that there be an even number that are 
signant—two indices in a 4-way or 5-way matrix, two or four indices in a 
6-way or 7-way matrix, and so on. We therefore lay down the following: 

Definition of a Determinant.—A determinant of a p-way matrix is the 
sum of all the terms that can be formed by taking a set of perjunctive elements 
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as factors and prefixing the product of the signs of an even number of chosen 
ranges. 

A permanent might be viewed as one extreme, where the even number is 
zero; but file-multiplication and dependent processes have no application to 
permanents, and so we prefer to mention them explicitly when a theorem is 
true with regard to them, and to understand that a determinant has at least 
two signant indices. 

By a full-sign determinant we shall mean a determinant as heretofore 
defined, with all, or all but one, of the indices signant, according as it is of 
even or odd class. 

If in any determinant two layers of a signant direction be interchanged, 
the sign of the determinant is changed. Hence, if two such layers are alike, 
the determinant vanishes. Hence a multiple of one such layer may be added 
to another without changing the value of the determinant. 

Both a determinant and a permanent are of course homogeneous linear 
functions of the elements of any layer and have the properties resulting. 


6. Decomposition.—A p-way n-layer determinant A can be decomposed 
into the algebraic sum of n! (py—1)-way determinants or permanents, as the 
case may be. Each of these components has n! components, and so on. UIti- 
mately we arrive at the expression of A as the algebraic sum of (n!)?—®? 2-way 
determinants or permanents. 

Let the matrix ||acg....,i|/@ of A be divided up into its (p—2)-way sub- 
layers of directions 1 and 2. Denote the sublayer common to the r-th layer of 
direction 1, and the s-th layer of direction 2 by d,,..5. ‘Take n perjunctive 


sublayers @qno....0) Var20....09 «++ +9 Gawno....0 to form a (p—1)-way matrix, 
Qa’10....0 
Qq''20....0 
Dainnd....0 


This is a component of the matrix of A; there are n! such components, and 
all of the locant of an element in A except the first index, is the locant of that 
" element in each of the (n—1) ! components in which it occurs. 

Denoting by the superscripts ~ and ~ the signant and nonsignant indices, 
and inserting a colon to isolate the locant of a component, we have: 


(1) oop. | P =U aoe... [P75 | 
B B 
a 


| (p—1). 
n ? 


(2) jacy:... |? =2Z+t. lag 
: (D,) 


ae 


(8) jaxg, [P=—=S. jaug.[ PP? 
a 





(4) lazy. [P=Z faaw.... 12. | 
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Each index beyond @ is signant or nonsignant on the right according as it is 
signant or nonsignant on the left. In (1) and (2), +, is the sign of the 
a-range in any transversal when the @-range reads 12... .n. 

Briefly, if « is signant (nonsignant) the components are signed (unsigned) 
and the signancy of @ is reversed (continued). 

We see that the components of a determinant may be determinants or may 
be permanents; but that the components of a permanent must be permanents. 

In verifying the formulas we must recall the fact that there are an even 
number of signant indices. Consider a term of |a,.7.,..|~ in formula (1). 
Let its elements be arranged so that the values of 6 are in the order 12... .n; 
then its sign is the product of the signs of the signant ranges beyond B. Pre- 
fixing +,, we find that we now have the sign proper to this term in |a@@,,,,|®. 
And this sign is preserved when the elements are permuted, provided that we 
make @ nonsignant, because there are an even number of signant ranges 
beyond ~. 

If the process be repeated to give the (p—2)-way components, we shall 
have (n!)*? matrices of the form 


a’ p'10....0 


Aq'1p1'20....0 





Dainpindnd. ee, 


To condense the corresponding formulas, we use a double superfix = =a 
the upper signs to be read together and the lower signs together in each 


formula: 


wa 


(1) | aB7....[2 = B tax 6] Gep.3... ie 
a, 





(2) \eag-.... 1° L+, toa. 1, 
e a, B (D,) 
(3) lazay....[P9= 2 xp|aep.5.... |e; 
a, 
(4) (a 1P=E laa... 
a,B j 


The formulas for complete decomposition, that is, separation into 2-way 


components, are: 
(1) (eV oti? =22 Ie... 201%": 
(2) ja... 291 =2 je...o018: 


(D,_2) 


+ a 
(3) la.ax [PHU (a. xx [2s 


+ + 
(4) la. ex[ PHU la... xx|2?. 


The “+” is the product of the signs of all the signant ranges before x. 
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It will be seen that the 2-way determinants in (1) and (2) and the 2-way 
permanents in (3) and (4) have for their rows files of A of the p-th direction 
(index A) i. e., rows of A; and that we get determinants when 4 is signant, 
and permanents when 4 is nonsignant. 

As we can put the indices of any determinant in any order before decom- 
posing it, the formulas are general with respect to such order. This remark 
applies to some of the later formulas. 











In general, if we wish to have r nonsignant indices a,a,....a,, and s 
signant indices 6,6,....@, come before the colon, an index y to come imme- 
diately after it (that is, to be the index of the range which is the base to which 
+, +2,)----, +,, relate), and wish ¢ nonsignant indices 6,6,....6, and wu 
signant indices ¢,«....¢, to follow y, the result is: 

WS 0880 A cece my tere eeee w | (r+st+t+u4+)D 
(eather Senay i=: Sey Oe n (D) 
aioe Ale coco meceen | (¢+44+) 
po x = pg, = p> oe. S me Way... 78y..0-Bs 27 Gaecee G6 €y 0000 Ey n ? 
gs ccOphieces+ Be 


where —" is to mean ~ or ~ according as wu is odd or even, regardless of 
whether y was originally signant or nonsignant. That is, the signancy of y 
on the right is to be so taken that there will be an even number of signant 
indices in the components. 

There is an interesting resemblance between the behavior of the signs 
~ and — and that of + and —. See (D,) and (D,), where, symbolically, 
Aas, AVE, SRE, Ve EY, A=, ete. And, generally, 
in (D), taking the superfixes of the a’s, the @’s and y, we find that 


— 


ane e me e 
wee we... AAS if s is even, and =~ if s is odd. 


7. Element-Multiplication—If from the elements of the matrices 


Whiriau..cmpllar’s || D6,65....B9ll 


we form a third matrix of class p-+-q—1 and order n, in whose rows the rows 
of the a-matrix and the rows of the b-matrix are found according to the pre- 
scription 


Ca,. . slips aPiye0+ Beat Fey... thy at By. «+ Bgl ? 


then it is seen that any transversal] of the c-matrix consists of a transversal of 
the a-matrix and a transversal of the b-matrix, every possible combination 
occurring just once. 

Let A be either a determinant or the permanent of the a-matrix, and let 
B be a determinant or the permanent of the b-matrix. Let C be a determinant 
or the permanent of the c-matrix, according to the result when the signancy of 
32 
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@,...-d, 1 in A, and of @,....@,_, in B is continued in C, and when u is made 
nonsignant in C if it is signant or is nonsignant in both A and B, but other- 
wise is made signant in C: 


~~ 5 _—_— ~~" 5 SONS TS OO net Se 


Then the evident theorem is: 
AR=C. 
The theorem includes as special cases both of the theorems of Section 3. 
In the case ~~=vL, if either A or B is of odd class, C has more than 
one nonsignant index. Thus determinants that are not full-sign determinants 
not only fit into the cases previously known, but also create a new case. 


8. Firle-Multiplication.—Given any two determinants with signant rows: 
AOR aaa BS OR aA | Ns 

let us compound every row of A into every row of B in the way familiar in 

the case of 2-way determinants and used in Section 4, so as to form a deter- 

minant C of class p+q—2 and order , according to the prescription 

re =2a,, poe eee ae 

that is, combine the locants of the rows of A and B to form the locant of an 

element of C, and continue the signancy of the indices in those locants. Then 


AB=C€. 


Proor: Completely decompose A and B into: 


ee sees | ee aor rs 
UE byes ss Eppa | OG....8.-2:BeaBel 5 
Da’, ....0’p211 Var, ....a’p212**** bz, silly tt Dey, veesBlq2 12° 
or, D+, Dar',....a! p21 Var's....ap 222 |? Lz Dgr,....A¢-221 D gy... .'%q-22 see 
Fm ete le sie ee: wibie tarares s eereee: hete hovers vere lee: asele evn sle (slen eco ete 


Multiply rowwise: 
daw, 000A po 1p Dey, oe Bl ag lh 2aw, 20008’ p-o 1p Dpr,....Bq-22u a ee 


AB=XL+,+5 Ldan,....a'y 2 be, — Ldarr,....ay 2 Dar, BM —-s2u tees 
A A 


Ca, eee Op» 1 Bi. . +B’ g-21 Cal, coe eA yo 1 B",.. - «Bg-o2 ee 


SS pS 
Bek, Catt... .a!y 22 B';...-B'q-ol Cathy. ...a!p 2 Py. P ggh * * * ° 





oeeeeeet vee e ee ee ieee eeeeeeeeeeseeeeeeevee 


The determinants in this sum are not components of C; but the mn! transver- 
sals of any one of them are m! of the transversals of C, and the (n!)?t**-n! 
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transversals of all are the (n!)?*?-* transversals of C. As to sign, a trans- 
versal 


, ; 4? 4? 
Ca’... 0, 91 BO... Br, Call... at! 280... Boy, 


has the sign +,+,-) +, in C, and the sign +, +, +, in this sum; and since 
there are an even number of signant indices in 6,8, ., the signs +,,-, and 
+, are alike. Therefore AB=C. 


rs Si wAA 3 —— YA AA 5 . 
ExamMPLe: A=|a i B= | byazaane|! ; 


0120 
CG ,a2B1b2838, Qa,az1 b6,6.6,8,1 7 Qa,a,2 b5,6,8,8.2 : 
(It will be convenient here and elsewhere to insert commas in- locants to 
bring out the structure of a matrix; they may be inserted, shifted or removed 
at pleasure, as they do not change the meaning of a locant.) 


os Ay, 11%, 12 Ae 1142, 12 sai Din, 11 bin, 12 Dio, 11 Diy0, 12 


— ’ — 


Ae, 21 F290 Ay, 9191, 29 Dove, 21 Doo0, 22 Dooy, 91 Boor, 22 


AB 1,1 Da, it+@y, 2 Din, 2 A411 Doo99, it 41,2 Doo99, 2 
As21 bin, 1+ Ae» Bun, 2 4221 Doo09, it Arg, » Doo00, 2 
41,1 D430, 1+ 41,2 D430, 2 A411 Doors, it Ay1,2 Door9, 2 


Qo9,1 By301, it Ag, 2 By301, 2 A211 Deore, 1+ Ag9, 2 Deore, 2 





C41, 1111 11, 2222 C41, 1121 C11, 2212 C 


Co2, 1111 C22, 2022 Coe, 1121 C22, 2212 





This proof depends on the rows of A and B being signant. Apart from 
the proof, it is plain that one of the two indices a,@, can not be signant and 
the other nonsignant, for that would give C an odd number of signant indices. 
And we proceed to show that both indices can not be nonsignant; from which 
it follows that a full-sign determinant will not serve to express the file- 
product of two determinants of odd class—the restriction stated by Cayley 
in announcing his law of multiplication. 

For, take n=2, to simplify the statement, In the transversal 


Ca),. ee +! pol Br. oe -B'g-2l Can,. o- A! 52 Bl", eee BM g-22 ? 
we find the monomial 
Dat,....alpgl2 Dey, veeBlqegl2 Vatt,....a!%y 022 Den, eB q222 9 


which does not consist of transversals of A and B, and we also find this 
monomial in the transversal 


Cy),. o* pol Bl,. oe -B!g-22 Can,. +. lly 2 Br. . ~Blg-gl } 





I 5 Es > AMES WIT Ng 


i 
4 
a 
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and these two transversals are affected with the same sign, if a, and 6, be 
assumed to be nonsignant, since there are then an even number of signant 
indices in a....a,., and in 6,....8,.,. Under this assumption, therefore, 
AB#C. The reasoning applies to every order of determinant and to every 
such monomial; if «, and @, in that monomial take one value v, times, and 
another value v, times, and so on, with v,+v,+....=m, then we shall find 
that monomial in v,!v,!.... transversals having the same sign prefixed. 

It is instructive to compare the matrices resulting from element- and file- 
multiplication. If the polynomial elements of the latter matrix be converted 
into nonsignant files by deleting the + signs, we have the former matrix, as is 
shown by the prescriptions. 


9. Crossed Decomposition.—Let us generalize, in form and in substance, 
the development of a determinant which is found in the proof of the law of 
file-multiplication in the previous section. Beginning with the decomposition 
of any determinant or permanent 

A=|dy,...04|2 
into 
Dt | Ga... ap ettp styl? 

in each component alike let us cause any set of indices before the colon, pro- 
vided that among them there are an even number of signant indices or else no 
signant indices, to take their nm sets of values not in the » rows but in the n 
columns (as in the c-determinants in the above proof). Let us call the new 
determinants or permanents, with the signs of the components from which 
they were formed, crossed components. Example, |az2,902,|$; a component, 
A411, 11 %11, 12 A411, 11 M122, 19 | : 
A222, 21 Ae22, 22 211, 21 Vo22, 22 

The sum of these crossed components is equal to A. For, let the trans- 








; the derived crossed component, 








posed indices be a,;....@)-2, and consider any transversal of A: 
Ga,. . «+ Cyuglty Oa, ...<tagu, °** Baf?.... af aag”. 
Let r,r,....17, be such a permutation of 12....n that 
ag? 1, ae 3, ... . «a =n. 


This transversal will be found once and only once, in and only in that crossed 
component whose main diagonal is 


Gay... 0,08, .... 08911 Gay’....a, af? ....08%22°>:* Bam... afPafte).... akan, 


i. €., in the crossed component 
Gay... 0, OP... Ql] Ga, ....a, af... .a8%12.... 


Gaz’... 05 af?,...-a0921 Gay... .a, aft... .08%22.... | 
Sao eae pic Cos HEMT aaadare Ghee ine tao eae see 

















ae 
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and the sign will be right, for if by the symbol 


i eS 
Qnr+1...:p—2 


we denote the product of the signs of the signant ranges among a,,;...-.@,~s 
when they take their values in the order indicated by the symbols (r,), (12), 
.., (7,), then we have the equation 


ul) tt Vidal. Ol el Sl ie =e Fas 
Generalizing in substance, let us decompose any determinant or perma- 
nent A into components of class q: 


= e ( 
A=xX+ | Qoy,...O1n, 421... Lop, » eee Ae Aoh,* Ayo. - + - Oh ©. 


q 


sees 


the indices before the colon being separated into q groups with an even number 
of signant indices (0 even) in each group other than the first. Alter each 
component by causing each group other than the first to take its » sets of 
values not in the » 1st-way layers, but in the m s-th-way layers for the s-th 
group. By this alteration the main diagonal is unchanged. Retain the signs 
of the original components for these crossed components; in other words, give 
each crossed component the sign of its main diagonal term as a term of A. 
The sum of these crossed components is equal to A. 


Exampie: The determinant | @¥0,3,2,2,0.aa4,|2 has 3-way components, 


3) 
= sya | 6 
x —— Gyg —— Agi —— Agg —— Ag, —— Age Wa, 104 202102003105, * G4deGg;2 ° 


Ae, 22, 22,, 211 Fee, 22, 22, 212 
As, 22, 22,, 221 Tee, 22, 22, 229 


441,13, 11,, 111 911, 11, 11,, 112 
41, 11, 11,, 121 %11, 11, 11,, 122 
ing crossed component: 


A component: + | ; the correspond- 








| 
Ae2, 11, 11,, 211 %22, 11, 22,, 212 
Geo 22, 11,, 221 V2, 22, 22,, 292 


441, 11, 11,, 111 411, 11, 22,, 112 
41, 22, 11,, 121 11, 22, 22,, 122 











n 
Proor: Let 1 dal? afPas? afP..ai? agp tas? af? be any transversal 
t=1 gi 


of A. Take g—1 permutations of 12... .n, Viz., rare... +m» SH=2,3,...-5 |g, 
such that 


af =1, ef? =2,. ...., 0a" ou 


The transversal will be found once and only once, in and only in that crossed 
component whose main diagonal is: 


n 
Ta (a?...afasie? af? aft? afsrét... .t). 
t= 

And the sign will be right, since, for each value of s, 


SEA s.5:8) SS Eis Pes 
Os1....sh, Aste ish 



























SL a marvin e baton are 


VEN kee hunt tei Gr eels Gab nf cele iae es 
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Thus, in the example above, the transversal 
112911121 V21129912 


must be found, if at all, in the position 


Ae, 11, 22,, 212 











yy, 22, 11,, 121 


in a crossed component; and there is evidently one and only one crossed com- 


ponent in which it is found. 
The development in crossed components may be written: 


Oy ++ + Ayn, Ay 

Gnsis)ss sie a 
A=S (tg. lad mem, % be 
( 8 ae | oe os wwe o> ee E 

Agus ++ +Agn, | Ag 


10. Raising and Lowering the Class: The “ Determinant-Permanent.”— 
One may arbitrarily raise the class of a given determinant or permanent by 
introducing one or more new indices whose values are determined by the 
values of one or more of the original indices, suitably adjusting their signancy. 
‘For example, using Kronecker’s symbol 6,..,=1 if i:=....=1,, otherwise 
=0, we have: 


Ay Ap | — | %1,1%1,2 


A21499 





legs P= lagesaa |? ; ie, 











Ao2,1 22,2 


And in general we may introduce as many nonsignant indices as we like, each 
having any one-to-one correspondence that we like with one of the original 
indices; and may then arbitrarily render signant any two of the entire set of 
indices which happen to have a one-to-one correspondence with each other 
such that both have the same sign. 

In particular, if a determinant of even class have two nonsignant indices, 
there may be introduced a signant index which shall everywhere take the same 
values as one of the nonsignant indices, the latter being also made signant, 
and thus being doubled, the result being a determinant of odd class with one 
nonsignant index. I. €e., 


lazy... 2°= [Saaaagy...1. 


It is by this kind of a determinant of odd class that the product of two full- 
sign determinants A® and B® by file-multiplication has heretofore been 
expressed, when p and gq were odd, as a determinant C of class p+q—1 (not 
p+q—2); the “fixed index” (nonsignant index) of A has been doubled in C, 




















eek ees 
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and the “ fixed index” of B has been made the “fixed index” of C. Of course 
C would consist largely of zeros.* 

On the other hand, let us start with a full-sign determinant A, of odd 
class, of a type of which C in the last paragraph is a particular case; wherein 
a group of indices take the same values, another group take the same values, 
and so on, there being r groups with an even number of indices in each, 
s groups with an odd number in each, and ¢ single indices, among the latter 
being the nonsignant index t,.f Give to the elements new locants by striking 
out all but one of each group of indices, put them into a matrix 


i >... ites eal oo 


and consider the determinant 


Mam a nts ~ | (r+s4+). 
——s Py cece PrOyoees TgTye0++ Tea Te | b 


noting that s+ is necessarily odd. Evidently, 
A'= A 


Now 4’ is not a full-sign determinant, but from A’ there could be formed 
a function not involving anything but full-sign determinants and permanents, 
viz., the function invented by Gegenbauer + and called by him a “ determinant- 
permanent.” Decompose 4’, using (D) of Section 6, into . 
(r+1) 


V+ ies ae ihisesl 
. bs 


fe | a; 
Og 8 8 8 Tig | My. 0 Og Tye eee Tens TePi- oes Py 





There are (n!)**’' of these permanents, and in them the (o,....¢,7,....%_4)- 
ranges are differently written before the t,-range, and prefixed to each of them 
is the sign which is the product of the signs of these ranges. Such being 
the case, we can arbitrarily construct a determinant B with purely formal 
elements: B= |bp0-957°3,5,/ 24° 


coe Sea 
which will have the property that if for each transversal in its expansion in 
terms we substitute the corresponding permanent, the resulting function will 
be equal to A. This function is a “ determinant-permanent” of class s+¢ and 
genus r+1. 

If A be made of even class by deleting t,, A’ will be of class r+s+#—l, 


its components will be of class r, the indices after the colon will be Q1----Prs and 
B will become 


lboa33 “s ore ’ 
Oye0ee Og Tye ees Te | NM ? 





*M. Lecat, “Sur la multiplication des déterminants,” Ann. Soc. Sci, de Bruwelles, Vol. XXXVII, 


Part 2, p. 285. 


¢ Abrégé, p. 32. 
+ L. Gegenbauer, “ Einige Siitze iiber Determinanten hiéheren Ranges,” Denkschr. Akad. Wien, Vol. 
LVII (1890), p. 735. Abrégé (see Note 3), p. 32. 
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that is, the “determinant-permanent” will be of class s+¢ and genus r. 

In either case the “determinant-permanent” is simply a decomposition 
of A’. 

It has been said that a “determinant-permanent” is necessarily of odd 
class.* But the decomposition formula (D) shows us that there: are two 
possible decompositions of a determinant |a3,..-a7,,,/:1:z,| into a sum of per- 


manents, such as to leave only signant indices before the colon, namely: 


ae coe =, Aq... .ay ry oe si ’ 
atte az, Gy.---Apzyi ay . wea 


Let A’, therefore, be decomposed into 


























sas a CHK), 
pe a i Ts a PS mee ey Tt ies 1(T¢)P4-- ye 
then 
a woe et 2. 
= =| o> LOeTy. hie 1 
and the “ determinant-permanent” is of even class. 
EXAMPLE: 
-— lananarnre oa I a ee B.S 
me lavage | {== (1) + 9 Gee 57 (P= (11) D+.|\0.: 557 |". 
as ntl 2 B= (ii) |baa|$ . 
(i) Replace }43,b20 + D432D001 — Dyno —Dy.0bo11 
+ +. + . F 
by Andie Gy3214 122 A211 i212 Aj2210 1209 
, , 
399103099 99119019 212142399 931149312 
+ > 
- a a ‘ 
that is, by 1, 111, 11,1011 4, 111, 11,2222} 4 ee eo 
Ae, 229, 22,1111 42, 292, 22, 2202 
(ii) Replace Dy,b50 — Dy bo, 
, , , , , , + 
by Q4111%1112. | 42211 F019 | __ | %111%2112 Geno. 
, , , , , 
312141192 | 42221 % 2299 A301 %o199 | @}291 44290 

















and translate into a’s as in (i). 
In the particular case of the determinant C described in this section, we 


have, first: 


| asa, LP Og aca P= | Cage e RRA 
=Utg. +++ ta, tp-+-- i) et 2 (1) 


which gives the product of A and B in the form of a “ determinant-permanent ” 
of class p+q—3 and genus 2, a form previously known. We attach the prime 





* Abrégé, p. 32. 
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to c to give notice that one index of the group a,a, has been dropped, just as 
a’ was used above when all but one index in each group had been deleted. 

And secondly we have the new form, obtained by writing after the colon 
any one of the indices that are before it in (1): 


tgs «i Mg tes 


2 





, vw ew (3) 
vee e Hg | Cas...dpBy-++-Bg-2!By-1Biay 1m 9 (2) 


a “determinant-permanent” of even class p+q—4 and genus 3. 


11. A Product Determinant.—A certain example given by Muir was 
generalized by Metzler in a paper “Ona Determinant Each of Whose Elements 
is the Product of k Factors.”* HE. H. Moore contributed to the subject + and 
mentioned that a particular case of the theorem was ascribed to Kronecker. 
All this work was in two dimensions. Then von Sterneck extended Kronecker’s 
result to p dimensions.t We shall extend Metzler’s theorem to p dimensions, 
including von Sterneck’s generalization as a special case. 

First let us deal with two sets of determinants, A, A®,...., A; 
BY, B®, ...., B™: 


AP x= | gOy sy er (P) . BO = |b@3°---59 ree 


see D7 Apat.+0+ Ay By..-- BgBgi1----Bg lt ° 


Form a p-way determinant A=AA®....A™, of order kl, by placing 
A®, A®,...., A® along the main diagonal, all other elements being zeros. 
Using the bipartite signs of order 


i eee oe: 2 ee ee eee 
such as Moore employs, we shall have the prescription 


Ga) (0B) 2278) (8 paaB) 22 (GB) = Oyo * 

In the same way, form a g-way determinant of order kl, with B®, 
B®,...., B™; then alter its form by a rearrangement of layers, placing the 1 
layers of the first direction containing B” in the positions 11, 12,...., 11, 
those containing B® in the positions 21, 22,...., 21, and so on, and do this 
for each direction. The resulting determinant B=B’B®....B® and the 
prescription is: im 

D a) 22GB) @Byasd 2 Bg) = OB. Bg" 





*Am. Math. Monthly, Vol. VII (1900), p. 151. 

+“ A Fundamental Remark Concerning Determinantal Notations with the Evaluation of an Impor- 
tant Determinant of Special Form,” Annals of Math., Vol. (2)1, p. 177. 

£R. D. von Sterneck, “ Ausdehnung eines Kronecker’schen Satzes auf Determinanten héheren 
Ranges,” Rend. Palermo, Vol. XXX (1910), p. 58. Lecat points out the fact that von Sterneck’s theorem 
does not hold when the classes are both odd; Abrégé. p. 63. The present extension is not subject to that 
restriction. 


33 
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Multiply together A and B by rows into a determinant U of class pt+q—2 
and order ki. The elements of U which are not zeros are monomials, since a 
row of A that is not blank contains nonzero elements only in the B-th set of k 
places, while a row of B that is not blank contains only one nonzero element 
in each set of k places. The a-element and b-element whose product is a 
u-element will be the a-element in whose locant a,=a and the b-element in 
whose locant 6,= 8; we therefore change a to a, and 6 to @, to form the 
resulting prescription: 


U (a,Bq)..- (GB) (A7asBy) 
where at least (a, ,6,) and (a,8,_,) are signant. 
This gives the theorem 
U=—=A” AVP BY Be 


For 2-way determinants the prescription becomes 


~~ ~ = (B2) (a, ) 
u (182) (4281) =a, bis, ’ 


which agrees with the theorem designated 7, by Moore. 
tA? es? —.... 2A? oad, coy, ee Bn PD? —.... =< BB, we 
have: 
Ux AP", 
which is an extension of von Sterneck’s theorem to less than full-sign deter- 
minants. 


ExampLe: p=3, g=3, k=2, 1=2. A’=|a'vaa,|?, A”=|a'vaa,|%, 





=< Ady 
B’=|b’xa3|%, B’=|b’xan|®. For A and B we have: 
(11) (21) : (12) (22) 

(11) (21) (12) (22) = (11) (21) (12) (22) = (11) (21) (12) (22) ~— (11) (21) (12) (22) 
(11) |] Qindne + + | Gndoe 
(21) | GierGiee + + | Aoe1Ao99 ; 
(12)p - se ep et et Poe + Gnas | + aanGee |’ 
(22), - 5 ts efor ttt ft + peice | + + Gee Ae9e 
(11) 111 * Die > My ou * Dore : 
(21) pi at, ‘te “o> * bin + Dore 
(12) Din ’ 129 +e oR Bly 292 ; 
(22)] > * + + | * Din + Duel * + + + ft Don + Date 

















The prescription is: 


~ ee eee eS (a3) 
u (83) (4283) (4381) (4382) — Qa,dyas DB:8.6, . 





























subdivisions. 


element under the a-factor. 
first index of wu appear in the four large horizontal subdivisions; those of the 
second in the large vertica] subdivisions; and those of the third and fourth in 
the horizontal and vertical lines in each square common to two intersecting 
Thus one of the sixteen rows of A is associated with each of 
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We shall condense the representation of U by writing the b-factor of each 


It is convenient to have the four values of the 


the sixteen squares, while one of the sixteen rows of B is associated with each 
of the sixteen places in each square (the same place in every square): 


(11) 


(21) 
(11) 
(12) 





(22) 


(11) 


| (21) 





(11) (21) (22) 
(11) (21) (12) (22) (11) (21) (12) (22) (11) (21) (12) (22) = (11) (21) (12) (22) 
, , , , 
Qin Qin Ay01 Ayo 
, / / / 
bin : Dior f Di : 121 
Ay12 Ar12 Aj29 Ay29 
4? bb”? a7 47 
111 121 111 121 
, , , , 
Qin Ay Ayo Ay21 
, , b; . , 
211 221 211 221 
Arie A119 Ai0 Ay09 
47 4? vr v7 
211 221 211 221 
, , , / 
A211 Qe Ago1 A291 
, , r b’ . . 
111 121 111 121 
, , , , 
As12 Aare A209 A029 
a? bi? v7 b”. 
J , wal 121 
Ae11 Qoy1 Aoo1 A201 
, , ’ , 
21 221 211 221 
A212 Ao12 Ao90 A202 
4? be ° a? 47 
211 221 211 221 
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We deal next with three sets of determinants: 


-. x, pas a. a, ee aad i neate a, wecuial a”, 

| amet _, errr oo; , Meets -, ee . ree 7, 

, , > aithaed Co, , prone Cc, eee of. eee as gm; 
where 


(a8) — rn: ~~wA (r) 
C lc > TeYaeas- at be 


Formms now U0? = AM? A™ .... AP BOE™ ....B™, and 20 01 up to 
y= A™...,. BB, construct C=U™....0™ (aa we did 4) ao that 


Tae Li -e ~~ sxe. ad — =a, » Og : 
(a, Bg) cons (a7¢Bq7) (Gp41B9) eves (ap-1Bq) (apBy'Y) ASse (apByY) (apBys1)-- @pBga) ~Bg ? 
and construct C (as we did B) so that 


— (af) 
aaa. ORE Oe 


© (ap y1).22. (GB Yn) @BYa43). 0s (@B Ye) 

Multiply together U and C by rows into a determinant V of class 
p+q+r—4 and order kim. The elements of V which are not zeros are 
monomials, and a,=a, 8,1=(@, y,=y. Detaching the suffixes of a, 8 and y in 


the locant of an element of V, we have the prescription: 


ws LF FLT pp pp ppp | 
OY Bs: 2.4 -++- 1....g g+1...g—2 q—1....¢g—1 q—1....q—1 
7t tf Cit on f..0f 1.8 h+1.... r—1 | 


=a (har ) » OB, citehes p. 


Here at least the indices (a, 4@,y,) and (a,@0,17,-1) are signant. The super- 
fixes on the right are the sets of two consecutive indices in the sequence 


B, Yr Boa : 
The prescription applicable to four sets of determinants is obtainable by: 
(i) changing v to w; (ii) subjoining to the locant the line 


oe ee See eee eee ee a 


the last s falling under y,_,, the 1 under y,_,, and the values p, gq—1, and r—1 
being continued over 6,....6,-;; (iii) annexing on the right 

AS40P ae) ‘ 
and (iv) inserting 4, after y, and before a, in the first three superfixes, so that 
the superfixes are now the sets of three consecutive indices in the sequence 


B, Yr é, Ap Bs Yr-1+ 


Hach prescription is obtained in like manner from the prescription that 
precedes, and we have the following general theorem: 
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The product of r sets of determinants 


AOI) == | gone Cre ceaaee, x Pass (Pr) 
as Any Ang-+++ Ant, nt +--+» Anp, | Mn 9 


where h=1, 2,...., r, and where g, 1s the h-th set of r—1 consecutive indices 
in the sequence 


A2po%3ps* + + + Hep, Ap, %2, py—1%8, pp—1> + + * Sr—1, pps 


is expressible as a determinant A of class py+....+p,—2(r—1) and order 
NyN,....2,, nm which all elements are zeros excepting those given by the pre- 
scription 


~ ~~ ee 


(G1191) ---- ( O1,, 91) (Otay. 43 oo (1,1 91) 


oeeeeeereeeeee eee eee ee eee eee ee eee 


-_~ 





L (1 Ir) + (np, Dr) ( Orp41Ir) +++ (Srp, 1 Ir) 
where k=2,3,...., r—1. 

Remark that at least the indices (a,,419:) and (a,),19,) are signant. 

If all the determinants are of two dimensions, the prescription takes the 
form 


“ 
a Hale 

(41 Ag0Ago. ++ Are) (Ay2M21g1.++-Gp_y53 Gry) Pr An34n2 ? 
~= 


the sequence to which g, applies being 
Ag9Ag0- “se K,9h42%o1 Hg - a -A, 41° 


This agrees with the theorem of Metzler designated 7’, by Moore. 


12. An Application to Transvectants.* 


TueEorEM: If the n?*? k-th transvectants of all possible pairs of binary 
forms taken from two sets 


—_— >. IIg,....6,116”5 n2k +2, 


be made the elements of a determinant 


Can lox-"4s a a a Cvceimy | (p+q) == odd, 

Oy e+ + Oy Agar ++ ApBy.-+- BrBrar----Ba n 9 q—h odd, 
= k 

with Ca,....ap Bee-+-Bg- Ma... Pz,. Ba) ; 

then C=z0. 





* Special cases of this theorem have been given by L. Gegenbauer (loc. cit.); see Abrégé, p. 94, 
where an important restriction on Gegenbauer’s results, not mentioned by him, is brought out. 





{ 
; 
4 
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Proor: We can form two null determinants A and B, whose elements are 
k-th derivatives of the /’s and ¢’s respectively, with suitable numerical fac- 
tors, such that 4-B=C. Letting m,,...0,, Mg,....6, be the degrees of fa,...4,, 


Pp,....6,) Tespectively, set up: 
(M,....:0,7 78) ! _ 


? Go,....0a,°™ 





_—_ were a eA +1 
A=|a Xr © ) 





Gyo eee Dy Ag41+0 + Ay csc ! = "i ? 
’ _—*.. k 
B=|by see sae | a+», i =(—(, k f (Ug,....6,—*) ! d Pp,....B, 
si 1 ++++PraPnii-+--Bq n LeeesPgn” iba i A—19k—A+1° 
[B,....Ba ! e oa 


These prescriptions fill only the first k+1 places in each row (file of the last 
direction), and we shall fill the remaining places in each row with zeros, the 
result being that one or more layers in A and in B will consist of zeros, 
whence A=0 and B=0. Obviously A+: B=C, the multiplication being of row 


into row. 
Corotuary 1. If the n* k-th transvectants of all possible pairs of binary 
forms taken from two sets 


| ees Pir- ++ Pry 


Tasso Tyg Pni- ++ Pan 
be made the elements of a 4-way determinant with two signant indices: 
C= | Ca,a2B; Be | “ws Caya28;82== (fara, ? 5,62) , 
then C=0. 
Corotuary 2. If the n’ k-th transvectants of all possible pairs of binary 
forms taken from two sets 
fis fey esssr tas Pry Par e+ +1 Pn» NZK+2, 


be made the elements of a 2-way determinant 


(frs-Pr)" + +++ (fry Pn)® 


Sas $i)". eee es ?,)* n 
then C=0. 


This corollary includes as a special case Gordan’s result: 


(fr i)" ++» (fas Oa)? 


(fs; 9)’. paises (fs, ps)” 4 
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On a Certain General Class of Functional Equations." 


By W. Haroutp WILson. 





§1. Introduction and General Considerations. 


Addition formulae of the general type 
G[f(a+y), f(x), f(y) ]1=0, 


where G is a polynomial in its three arguments, play a prominent role in the 
theory of elliptic functions. A natural generalization of such formulae is 


P[a, y, f("), f(y), f(aet+Biy), ---+) f(ee+Bny) ]=0, (I) 


where 

(i) P denotes a polynomial in its n+4 arguments such that every argu- 
ment involving f is explicitly present ; 

(ii) # and y are independent variables; 

(iii) a, and 8;, i=1, 2,...., , are given constants;+ and, 

(iv) f(x) is an unknown single-valued function to be determined so that 
equation (I) shall be identically satisfied. + 


Equation (1) is said to be of order n. The degree m of P in the function 
f is said to be the degree of equation (1). 


‘ 





* Read before the American Mathematical Society (at Chicago), April 6, 1917. 

+ For the purposes of this paper it is convenient to carry certain hypotheses in regard to the a’s 
and f’s. A statement of these hypotheses is to be found below. 

{A theorem of some interest in the general theory of these functional equations is that every solu- 
tion f(x) of equation (1) is a solution of a similar equation in which x and y oceur only in the arguments 
of the function f. To prove this arrange P as a polynomial in # and y. The substitutions e=s-+ k,t, 
yt, where k,=0 and 

— rn  » Ban _ Bs 

kpk+t =" kik t+ <a ’ 
A=0,1,....,1—1, h, j=1, 2,...., , transform (I) into equations similar to (I) such that the highest 
degree in s and ¢ is the same for all of them. It is easily seen that a finite number of non-zero k’s may 
be employed such that the variables s and ¢t may be eliminated from these equations, in so far as they occur 
as coefficients, by Sylvester’s dialytic method of elimination. The result of this elimination is an equa- 
tion (II) which states that a polynomial Q in the function f has the value zero. The arguments of f are 
linear combinations of two independent variables s and ¢. Hence (II) is similar to (I), although in 
general its order and degree will differ from those of (I). Furthermore, if no two arguments of f in (I) 
are proportional, then no two arguments of f in (II) are proportional. 
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Cauchy * discussed two special cases of (I) and two related equations, 


— f(ety) =f) thy), fa+y)=f(a)f(y), 
f(ay)=f(x)+f(y),  f(xy)=f(x) f(y). 

One or the other of the first two of these equations has since been treated t by 

Darboux, E. B. Wilson, Vallée Poussin, Schimmack and Hamel. Carmichael { 

has given a generalization of the Cauchy equations while Jensen§ has dis- 

cussed several applications of them. Cauchy || has treated the equation 


$(x+Yy) +$(x—y) =29 (x) p(y). 
Carmichael ff has considered the equations 
h(a+y)h(a—y) =h* (x) +h*(y)—e’,  g(at+y)g(a—y) =9' (x) —9?(y). 
Van Vleck and H’Doubler** have discussed the equation 


v(at+y)¥(z—y) =[d(a)d(y) ]?. 
Other related equations have been considered by several writers, and systems 
of functional equations have also been treated. 

It seems that no systematic account of a general theory for equations of 
the form (I) has ever been undertaken. This paper is designed to contribute 
to such an account. The equations considered in the principal part ($§ 2 to 8) 
of the paper are linear homogeneous equations with constant coefficients, 
They may be written in the form 


LD yst (ae +Biy) +7 usaf (@) +Yns2f (y) =0. (1) 


It will be shown that if some a’s and $’s having different subscripts are zero 
and no ratio a;/@; of non-zero a’s-and #’s is distinct from all the remaining 
ratios, the equation is exceptional. The exceptional case receives mention 
only in§$6and11. There is no loss of generality in assuming that no a is 
zero in the non-exceptional case. For convenience in exposition the hypothesis 
will be carried in the text that in addition to no « being zero, no # is zero, and 
no two ratios a;/8; are equal. The additional argumentation for the remaining 
non-exceptional cases is supplied in footnotes. 





* Cours d’Analyse (1821), Chapter 5. Cauchy treated the last two equations by transforming them 
into the first two equations. It is obvious that similar transformations may be applied to reduce more 
general equations to the form of those considered in this paper. 

+ Darboux, Mathematische Annalen, Vol. XVII (1880), p. 56. E.B. Wilson, Annals of Mathematics, 
Vol. I, Ser. 2 (1899), p. 47. Vallée Poussin, Cours @Analyse infinitésimale (1903), p. 30. Schimmack, 
Nova Acta, Vol. XC, p.5. Hamel, Mathematische Annalen, Vol. LX (1905), p. 459. 

{ American Mathematical Monthly, Vol. XVIII (1911), p. 198. 

§ Tidsskrift for Mathematik, Vol. II, Ser. 4 (1878), p. 149. 

|| Cours d@’ Analyse (1821), p. 114. 

{ American Mathematical Monthly, Vol. XVI (1909), p. 180. 

** Transactions of the American Mathematical Society, Vol. XVII (1916), p. 9. 
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A normal equation of order m is derived (in §2) which is satisfied by 
every solution of any non-exceptional equation (1) of order n. This normal 
equation forms a foundation upon which the entire development of the theory 
of equation (1) is based. Any normal] solution may be uniquely determined 
at each point of a dense set covering the complex plane if it is given at the 
vertices of a certain triangular network (§§3,4). It is shown in $5 that the 
normal solution analytic in the neighborhood of the point zero of the complex 
plane is an arbitrary polynomial in x of degree n. It is also shown that the 
normal solution continuous im the neighborhood of the point zero of the com- 
plex plane is an arbitrary polynomial in u and v of degree n where u and v are 
real and <=u+vV—1. The normal solution analytic along any line in the 
finite complex plane is also an arbitrary polynomial in x of degree n and the 
normal solution continuous along any line in the finite complex plane is an 
arbitrary polynomial in wu of degree n if the line is not parallel to the axis of 
imaginaries and an arbitrary polynomial in v of degree n if the line is not 
parallel to: the axis of reals. The analytic and continuous solutions of (1) are 
found ($6) from the normal] solution. Examples are exhibited in $6 which 
show that equations of type (1) may have non-trivial continuous solutions, but 
no non-trivial analytic solutions, while other examples show that equations of 
type (1) may have analytic solutions which are also the most general continu- 
ous solutions. A converse theorem is briefly considered in ¢7. It is shown 
(in $8) that if a function f(x) satisfying an equation of type (1) has a point 
of discontinuity in the finite complex plane [or on any line in the finite complex 
plane] it has a point of discontinuity in every region [interval] of the plane 
[line], however small. 

Equation (1) is employed ($9) to solve certain equations of the type 


Le(, y)f(a¢%+ By) HOnir (2, y) Ff (4) +O,42(%, YF (Y) + On43(%, y) =, 


where the @’s are known functions. Equation (1) is also employed (§ 10) to 
find all analytic solutions, and in some cases, the continuous solutions of 
binomial] equations of the type 


1 [f(a,7+6,y) = Th [f(ae+By) |* [f(y) 1%, 


where no a is zero, C is a constant and the y’s are constants of which the real 
parts are positive. Pexider* used the first Cauchy equation to solve 


f(x) +o(y)=¥(a+y). 
In $11 the method of obtaining the solutions of (1) is used to solve the 
equation 


M= 


vif (ae7+B,y) +Ynsifnss (2%) +Ynsofnse(Y) =0, 


t=1 





* Monatshefte fiir Mathematik und Physik, Vol. XIV (1903), p. 293. 
34 
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of which the equation considered by Pexider is a special case. It is proved 
that when no two arguments of f/ in the foregoing equation are proportional, 
each continuous solution f is a polynomial of degree not greater than n. 


§2. Reduction to a Normal Form. 


The solution f(x) of the general n-th order equation 


Dif (ae+ Bey) +Y neal (®) +Ynsaf (y) =0, (1) 
is contained in that of an equation of the same form and same order [equation 
(6) below] in which each a, 8, y is a given integer. The derivation of (6) 
from (1) is accomplished by elimination. 

If (1) is subtracted from the equation derived from (1) by replacing y 
by y+#,41, the result is 


Zylf (e+ Biyt Biter) —f (a+ BY) +7 ns2lf (y+ ter) f(y) 1=0. (2) 
If (2) is subtracted from the equation derived from (2) by replacing x by 
x—B,t, and y by y+a,t,, the result is 

Lyilf (wet Biy + Anti + Bitner) —f(a¢+B,y+ Arts) 


—f(aet+Byt+Bitra) +f(aet+By)] 
+Ynsolf(ytoyts+ti41) —f(y+a,t) —f(yt+tass) +f(y) }=0, (3) 


A= 4,8;—4,8;. 
It is easily seen that this is true because the given substitutions leave a,x%+ iy 
unchanged, but replace a,v+ 6;y, 11, by 
at+B;y+o,8t,—a:Pih=act+By+tAyt. 

In general, if an equation (a) results after such eliminations, then each 
argument of f in (a) is a linear expression in 2, y and certain ¢’s, the sub- 
scripts of the ¢t’s corresponding to those in the terms eliminated. The substi- 
tution of 


where 


a—6,t;forz and y+a,t, for y (4) 


gives rise to an equation (b) which differs from (a) by having each a,7+6;y 
of (a) replaced by a4,2+6;y+A,t;. Since A,=0 and since (4) does not affect 
the ¢’s that are found in (a), it follows that the difference formed by sub- 
tracting (a) from (b) contains no term for which i is equal to the fixed 
integer 7. Since n is finite these eliminations may be continued until only 
terms having the coefficients +y,,., remain. Moreover, the order of elimina- 
tion of the terms for which i=1, 2,...., n, is immaterial. 
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It is obvious that the equation resulting from these eliminations is linear 
and homogeneous. Since y,,.0 by the assumption that (1) is of order n, 
the equation may be simplified by dividing by y,,.. It is easily seen that 
there are (n+1)!/k!(n+1—k)! distinct terms in which the coefficient of f is 
(—1)*, and in which the argument of f is obtained by omitting k terms after 
the first from y+a,t,+at.+....+a,t,+¢,,,;. Moreover, no other such terms 
are possible. This is true for k=0,1,....,2+1. The substitutions ¢, for 
a,t,, t, for agt,,...., t, for a,t,, tar. for y, serve to completely determine an 
equation independent of the original a’s, @’s and y’s. If X, denotes the sum 
of the (n+1)!/k!(n+1—k)! terms in which the arguments are formed by 
omitting & t’s from X7*}t,, then the equation which is satisfied by f(z) is 


t=1 


Yo + Le - oe ~+ (—1)"%, + (—1)*'d,.1=0. (5) 
Equation (5) involves »+2 independent variables. In order to obtain a nor- 
malized equation having the same form and order as (1), let 4} =f,=....=t41 


=a and t,,,.=y, whence (5) becomes * 
: ! 
fi(mt+laty] +... + (aye EE pe ty) 
+....+(—1)"f(a@+y) + (—1)***F(y)=0. (6) 

From the foregoing considerations we see that every solution f(x) of equation 
(1) is a solution of equation (5) and of the normal equation (6). 

While the solutions f(a) of (1) are included among those of (6), it is not 
necessarily true that the solutions of (6) are included in those of (1). The 
following example suffices to prove this statement. Equation (6) for n=2 is 


f(3x+y)—3f(2a+y) +3f(x+y)—fly) =0. (7) 
As will be shown in § 5, the most general continuous solution of (7) over the 
finite complex z-plane is 


f(a) =aw+buv+ev?+dutevtf, . (8) 


where rx=u+vV—l1 (u and v real) and a, b, c, d,e and f are arbitrary con- 
stants. From the above considerations we see that the solution of any second 
order equation of form (1) is included among those of (8). However, substi- 
tution shows that for the equation 


f(2u-+y)—2f (w@+y) —2f (@) +f(y) =0 


* The results in (5) and (6) can be rendered more precise in special cases. If B;/a;—=Bn/an, Ajr=O 
and the substitution of «—,t,; for « and y+ a,t; for y leaves a,#+ Bry, as well as ajv+ B,y, unchanged. 
Hence the elimination of terms for which i=j also eliminates all terms for which ih where h is any 
value for which A,;,—=0. But when Aj;,=0, B;/a;—=£n/an, and therefore this proportionality is a neces- 
sary and sufficient condition for the simultaneous elimination of terms for more than one value of i. 
Hence the solution of any m-th order equation (1) in which no a is zero, is contained in the solution of 
equation (5) or of equation (6) where n is the number of distinct ratios B;/a; in (1). 




















268 Wiuson: On a Certain General Class of Functional Equations. 


we have d=e=f=0. In fact, such an equation need have no other than the 
trivial solution f(2)=0 as is shown by the equation 


f(2a+y)—f (x+y) —2f (xa) +f (y) =0. 


§3. Determination of the Normal Solution at the Vertices of a Network. 


The function f(x) to be considered in this section and in $$4, 5 is the 
solution of the normal equation developed in $2. It is to be observed that if 
f is given for the arguments ky, m+ky, 2um+ky, ....,uX%+ky, then f is 
known for the argument (n+1)x%+ky, by (6). By putting y=ju,+ky, for 


successive values 1,2, ...., of 7, f is determined by (6) for any argument 
mxz,+ky,, where m is any positive integer. Similarly, by successively giving 
to j the values —1, —2,...., f is determined by (6) for any argument 


mx,+ky,, where m is a negative integer. By interchanging 2 and y in (6) it 
is easily seen that f is known for all arguments ma,+ky,, where k is any 
integer or zero, if it is given for the arguments mx,, m%+Yyo, M%+2y,...., 
mx +ny,. Hence, if f is given for the arguments ma,+ky,, m,k=0,1,....,%, 
then it is known by (6) for all arguments ma,+ky,, where m is any integer 
whatever and k=0,1, ....,”, and finally, it is known by (6) for all arguments 
m2z,+ky,, where m and k are any integers whatever. 

It will now be proved that if f is known for the arguments ma+hky, 
where m and k are zero or positive integers and m+k<n-+1, then it is also 
known for the arguments ma,+ky,, where m, k=0,1,....,”. The figure 
illustrates the case for which n=4. The points designated by the small circles 


Immaginaries 











represent the arguments for which f is supposed known. The parallelogram 
formed by the lines joining the points 0, nz, na+ny,, and ny, will be denoted 
by x. The parallelogram x contains, on and within its boundary, all the points 
MxZ+ky, m, k=0,1,...., , and no other points as vertices. It is required, 
therefore, to find f at the vertices of ~ not designated by small circles. 
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Equation (5) furnishes a simple means of solution of the present problem. 
Let t,t. ...=t, = %, toys=....=t y.=—% and t,,,=0. By these substi 
tutions every argument in (5) is of the form mx,+ky, where m and k are 
zero or positive integers and m+k<n-+1, except in the first term in which 
m+k=n+1. Moreover, m<q and k<n+1—q; hence f is determined for the 
argument g%+ (n+1—q)y in terms of linear combinations of its values for 
arguments represented by vertices within or on the boundary of the parallelo- 
gram bounded by lines joining the four points 0, qx, gqx+(n+1—q)H, 
(n+1—q)y, (represented in the figure by dotted lines for g=2). By giving 
q the values 1, 2, ...., m, in succession, f is determined at each vertex of a 
which lies on the line joining the points 2%+ny) and na-+y, (the dot-and-dash 
line in the figure). It is to be noticed that each determination is made from a 
linear equation in one unknown with unit coefficient. Hence, the functional 
values so found are unique. 

Having determined f for the arguments ma,+ky, in a for which m+k 
=n-+1, it is easy to determine f at the vertices for which m+k=n+2. In 
(5) let t:=t,=....=t,=m% and t,,:=....=t,.2s=y%. Thus f is determined 
for the argument gz,+(n+2—q)y in terms of arguments represented by 
vertices within and on the boundary of the parallelogram formed by the four 
lines joining the points 0, qa, g%+(n+2—q)y%, (n+2—q)y. This deter- 
mination is also made by means of a linear equation involving but one 
unknown with the coefficient +1. Hence, no indetermination can be intro- 
duced. By giving g the successive values 2, 3, ...., , f is determined for all 
arguments in ~ which lie on the line joining the points 2a,+ny, and n%+2y. 

Proceeding in this manner f is determined for all the remaining argu- 
ments of 2 by giving t,,, the values 2y,, 3y),...., (n—1)y successively. 
For each value hy, of t,,., f must be determined for all arguments 

Quy+ (nth+1—q)y, gq=ht+l1, h+2,....,%, 
before giving t,,, the value (h+1)y,. As before, each determination is made 
by a linear equation in one unknown with coefficient +1. Hence, f is uniquely 
determined at all the vertices of 2, and we have the following result: 

Every solution f(x) of the normal equation (6) is known for the points 
ma2,+ky, where m and k are any integers of zero if it is given at the points 
mx+ky for which m and k are positive integers or zero and m+k<n-+l1. 


§4. Determination of the Normal Solution at a Dense Set of Points. 


It will now be shown that if the solution of the normal equation is known 
at all points mz,+ky,, where m and k are integers, then it may be found for 
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4mxa+4ky, and finally for 2~*mx,+2—‘ky for all integers m,k and s. If f is 
found for certain appropriate linear combinations of 2~*z and 2~*y,, then it 
is evident that 2-‘x, and 2~*y, may be regarded as were x and y,in§3. Hence, 
if it is most convenient to determine f at 2~*ma,+2~*ky,, m, k=0,1,...., 2m, 
it is also sufficient, in view of the argument of $3, for f will then be known at 
these points for all integers m and k. Furthermore, if it is shown that f is 
known at 2~'ma,+2—ky, for all integers m and k when it is known at ma+ky, 
then by regarding 2'~‘z, and 2'~*y, as % and y for the successive values 
s=2,3,...., it is seen that f is known at 2~*mxz,+2“‘ky, for all integers m, k 
and s. Therefore, since the points 2~*ma+2~‘ky,, m, k and s any integers 
whatever, form a dense set, it is sufficient to show that f is known at 
4ma,+tky,, m, k=0,1, ...., 2n, to show that f is known at all points of a 
dense set. 

Supposing f known at 4ma+sky, m, k=0, 2, ....,2n, it is only required 
to learn f at 4ma,+tky, m, k=1, 3, ....,2n—1. Replacing w by 3% and 


‘'y by that+ky, in (6), » equations are obtained by giving h the values 


n—1,....,1,0. The equations are linear inf. Alternate terms involve the 
arguments 4ma,+ky,, m=1, 3, ...., 2n—1, for which f is unknown. The n 
equations may be considered as linear equations in » unknowns. It is easily 
seen that the determinant of the coefficients of the unknowns is (—1)’A,, 
where v is n/2 or (n+1)/2 according as n is even or odd, and A, is a deter- 
minant of binomial coefficients such that the element in the i-th row and j-th 


aoe (n+1) !/(2j—i) \(n—2j +641) |, 
unless 27<2 or 27>n+i+1, in which case it is zero. Subtracting the 
(i+1)-th row from the i-th row, 7 having the values 1, 2, ...., »—1, in order, 


the element in the i-th row and j-th column, in, becomes 
(n—4j+2i1+2) (n+1) !/(27—1) !(n—2j +142) !, 
unless 2/<i or 2/>n+i+2, in which case it is zero. Adding to the j-th 


column the sum of the preceding columns, giving 7 the values n, n—1, ...., 2 
in order, A, assumes a form in which the element in the 7-th row and j-th 


column, in, is 


S=E (n—4h4 242) (m+1) 1/(2h—i) ! (n—2h+i+2) 1. 
R==1 


If for any given value of j this sum is n!/(2j--7) !(n—27+7)!, unless 27<1 
or 2j>n-+i, in which case it is zero, it is easily shown that for 7+1 it is 
n!/(2j7—i1+2) !(n—2j7+i—2)!, unless 27<i—2 or 27 >n+i—2, in which case 
it is zero. For 2j—i=0 or 2j;—i=1 the value of S is merely the first non- 
vanishing term, that is, 1 or , respectively. Since n!/(27—7) !(n—2j-+17) ! is 
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equal to 1 for 2;—71=0, and equal to n for 2;—i=1, it follows by induction 
that S has the value G=(n)!/(27—i) 1(n—25+4) 1, 


unless 27<i or 27>n-+471, in which case it is zero. For the n-th column the 
value of S is always zero for in since 27>n-+i%. For i=n, it is clear that 
the elements are affected only by the process of addition, and since every 
alternate term of (1+1)"*’ is involved, the last term in the n-th column is 2". 
Hence, A, has been so transformed that the principal (n—1)-rowed minor 
found by deleting the last row and last column of A, is A,_,, and the last 
column consists exclusively of zeros with the single exception of the element 
2” in the n-th row. Therefore A,=2"A,,. It is evident that A,—2'*? (A, is 
trivially 2). Therefore A,=2'*?*:-+"=2#"@+) which is distinct from zero for 
all finite values of ». Since (—1)’A,+0, it follows immediately that f is 





uniquely determined at the points $ma,+ky,, m=1, 3, ...., 2n—1, and hence 
for wai, 2, ...., 2m. 
Having determined f at the points 4ma,+ky,, m=1, 2,...., 2n, it is only 





necessary to let v=4y, and y=thy,+4mz2, in (6), letting h take the values 
n—I1,....,1, 0, to determine f at the points $ma+4ky,, m, k=0,1,...., 2n, 
for (—1)’A, is again the determinant of the coefficients of the unknown terms. 

Combining the results of §§3 and 4, it may be stated that f is determined 
/ at each point of the dense set 2-*mx,+2~“*ky, tf it is known at the points 
| mx%+ky, for which m and k are positive integers or zero and m+k<n-+1. 


§5. Solutions of the Normal Equation. 
Equation (6) may be written in the form 
att (n+1)! 
x (—1)’- 
an’) (nF) 
where a=u+vV—1, y=s+tV—l, and u, s, v, t are real. Let us seek a solu- 


tion f(x) of this equation in the form 
N h 
f(a) = 2B Degurors (9) 





fl (ju+s) +V—1(jv+t)]=0, 


in which the coefficients c,, are constants. Putting this value of f(a) in (6) 
we have 


ng . (n+1)! ae ’ 
Ee, Fe 00d hone , q( 47 h—q—(). 
2 ( 1) q!(n—j+1) ! 2 Bea (iu+s) (jvu+t) 0 
The terms involving u’~?s?v""*—"t" are 
n+1 . (n+1)! je ths q! (h—q)! —P eP),h—-q— 
ay h—p—r : I—P oP, q—Ter 
S| ecg? | pga aca 











j=0 

















272 Wuiuson: On a Certain General Class of Functional Equations. 


Since no term outside of the brackets involves j, it is evident that the given 
expression vanishes, for non-zero c,,, u, 8s, v and ¢, when and only when 
(n+1)! jo 
gi(n—j +1)! 
vanishes. Except for sign the quantity B,, is the (n+1)-th difference of 
ag’? for ~=0. The degree of each difference is one less than that of the 
preceding difference. Therefore B,, is zero when h—p—r<n. The particular 
value h obtained from h—p—r by putting p=r=0 must be included in the 
discussion of B,,. Since N was taken as the largest value of h, it follows that 
f(x), as given by (9), satisfies (6) if N=mn. Moreover, the value of c,, is 
arbitrary. 

Equation (9) shows that there are X7_.(h+1) =4(n+1) (n+2) arbitrary 
constants c,, which may be assigned at will. By §$3 and 4 it has been shown 
that when f is known at the 14+2+....+(n+1)=4$(n+1) (n+2) points 
ma+ky,, m and k positive integers or zero and m+k<n-+1, then it is known 
over a dense set of points covering the entire finite plane provided x and y 
are not collinear with the point zero. Since f(a), as given by (9), is continu- 
ous, it is only necessary to prove that each c,, is uniquely determined by 
assigning f at the given points mz,+ky, to know that f(x) is the most general 
continuous solution of (6) over the finite complex plane. This can be done by 
direct substitution,* but inasmuch as the determinant so formed is unwieldy, 
it is more easily accomplished by observing that the properties sought for any 
desired oblique network are readily deduced by a projective transformation 
from similar properties of a square array, on the axes of reals and imaginaries 
with the units 1 and 7, provided only that 2 and y, do not lie on the same 
straight line through the zero-point. Confining attention to the rectangular 





n+1 
bu © [ (8) 
j=0 


array mentioned, write a 

f(r) =X TAgu?v*” 

h=0 q=0 
where 
u?=u(u—1)....(u—q+1). 

For any integral value of u less than q,uw”=0. Beginning with =0 and 
proceeding outward through a triangular network similar to that employed in 
§3, it is possible to determine an A,, with each point of the net. Having 
determined the A,,’s, it is only necessary to expand and collect the terms of 
the expression for {(#) and compare with the expression involving the c,,’s to 


completely determine each c,,. 





*The value of the determinant of the coefficients may be shown to be 


n—1 
(VoSo— Uoto) M(M+D (n+2) TT [ (n—a) 1]2(a+1), 
a=0 
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Thus we see that the most general solution f(x)=f(u-+iv) of equation 
(6) continuous over the finite complex x-plane is an arbitrary polynomial in u 
and v of degree n. 

The analytic solution of (6) over the finite complex plane is that special 
case of the general continuous solution for which 

f(v)=a+aa+....+a,a*+..... 
Replacing x by u+iv, it is seen at once that a,u* must be zero for k>n, and 
hence a,=0,k>mn. Moreover, a,=¢yo, k<n. Hence the most general analytic 
solution of (6) is an arbitrary polynomial in x of degree n.* 

To obtain the most general continuous solution of (6) along any line in 
the finite complex plane, it is only necessary to observe that by the argument 
of §$3 and 4 it was proved that f is known at a dense set of points on the line 
if it is known at n+1 points of the line which are separated by some con- 
venient unit. The argument at the beginning of this section shows that along 
any line not parallel to the axis of imaginaries, 


f(x) es 


where each a is arbitrary, satisfies (6). Since f(x) is continuous, it remains 
only to show that the a’s are determined by the functional values at the »+1 
points on the line to know that f(x) is the most general continuous solution of 
(6) along the line. Substitution shows immediately that the a’s are uniquely 
determined by the n+1 values of f on the line. Hence the most general solu- 
tion of (6) continuous along any line not parallel to the axis of tmaginaries is 
an arbitrary polynomial in u of degree n. 

Similarly, the most general solution of (6) continuous along any line not 
parallel to the axis of reals is an arbitrary polynomial in v of degree n. 


au’, 


ita 


§6. Solutions of the Original Equation. 


The determination of the existence of any solution f(z) of the equation 
Lyf (ae+Biy) +Y ef (®) +Yus2f (y) =0, (1) 


and the determination of f(z) if it exists, is accomplished by substituting the 
corresponding solution of the normal equation (6) in (1). It has been shown 
in $5 that the general solution of (6) analytic over the finite complex x-plane 





* That the analytic solution f(),if it exists, is a polynomial of degree not greater than 2, is readily 
seen by direct differentiation of (1). If differentiation is made with respect to each of the variables y, 
ajv + By, j—1, 2,...., it results that f(*)(y)—=0 whenever k>n, because the arguments are inde- 
pendent in pairs which are not proportional. 


35 














274 Wruson: -On a Certain General Class of Functional Equations. 


or along any line in the finite complex x-plane is an arbitrary polynomial in x 
of degree n. Hence, the determination of the general solution of (1) analytic 
over the finite complex 2-plane or along a line in the finite complex 2-plane, if 
it exists, is accomplished by substituting 


f(@%)=aqtea+....+¢,2" 


in (1). Since the result of this substitution is an identity, the coefficient of 
any power of the variables must vanish. 

It is convenient to consider the terms of degree nm independently of the 
remaining terms. These terms are given by 


Ch [Dy (aetBy)"+ yur +7000") =0. 
Now ¢, is necessarily zero unless 
[Dy (a0+B:y)"+Yn410"+Yn42y"] =0. 
Placing the coefficients of this identity equal to zero, we have 
Baty yasr=0, py ai "Bry. =0,  Btyt+Ya2=0, k=1,2,....,;0—1. (10) 


Since only the ratios of the y’s are significant, it implies no loss of generality 
to assume y,,.—=—1. Under this assumption equations (10) may be employed 
to express the remaining y’s in terms of the a’s and @’s, provided f(x) con- 
tains a non-vanishing term of degree n. If we write r;=a;/@;, then for 


i<n+1 


n 
YiM1"2..- 1, / Bir (r,—1;), 


where the prime indicates that h does not take the value 7. Solution also gives 
Yn+1— (—1)"ryre. $6 Tae 

Therefore, if c,x", ¢,-0, is a term of the analytic solution f(x) of (1), it is 

necessary and sufficient that (1) may be written in the form 


2 n f(aat By) +(—1)"yre....7, f(x) —f(y) =90, 
Bir (n—r) 
h=1 
where 1;=4a,/8;. 
If the solution of (1) includes the term c,,2”, m<n and c,,0, then 


toss sef, 





Vary. t7au=0, : a”—*3ty,=0, Bry +Yn42=0, k=1,2,....,m—1. 
t=1 t=1 ¢=1 


Since there are but m-+1 linear equations in the y’s, they may be used to 
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express m-+1 y’s in terms of the a’s, 6’s and remaining y’s. For the first 
m-+1 y’s in terms of the remaining quantities, solution gives 

n m+1 
ae 2 BR (r;—T,) i. (—1)™* rr. eee Vea" j41- eee Vm+1Y n+2 


Y= : m+1 » (11) 


br IT’ (r;—1;) 
h=1 


where the prime indicates that h does not take the value 7. Hence, if c,,x7”, 
m<n and c,, 0, is a term of the analytic solution f(x) of (1), it is necessary 
and sufficient that each of the first m+1 y’s has the value given in (11). 

The computation involved in finding the most general solution of (1) con- 
tinuous over the entire finite plane is so tedious as to make it expedient to give 
results only for the general second order equation. The problem for any 
equation is merely a matter of substitution and algebraic computation. For 
the second order equation 


mill (yy + Qygt) B+ (04, + Dt) y] 

+f [ (oy + Aaqt) B+ (Ba; + Dot) y] +ysf (x) +ysf(y) =0, 
where 443, Qy25 03, B42, Ge, Gee, 021, Dog are real and i= V—1, the normal solu- 
tion with which substitution must be made is 

Cy -F Cypt + Cy 0 + Cag? + Cy UV + Coed”. 
This substitution shows that c, may be assigned different from zero when and 
and only when 7,+72+y7s+y.=0. It also shows that if we write 
Ay=Y10y +724 +s, By=Y1dy2 + 72002 Cy=Yibutyobatys, D,=y1b.+ Yebe2, 
then for c,, and ¢, to be independently arbitrary A,=B,=C,=D,=0. How- 
ever, Cy=Cyk, if A,=B,k,, Cy=D,k,, B, or D0, and ki=—1. - Further- 
more, if we write, 





A =y,0}, +720 t+Ys, B =A Ay, + Y2A01A22 y C=y Vig + 720% » 

D =yibi,+ybatys 1 E =Y10y b+ Y2Do1D09 , F=y7,bi,.+72b% ’ 

G=710 1). +7 201)0 , H=yAyby+y2dabe ’ J=Y10 2b. +7 2Aeeb00 5 

K =A 2b, +7 2022bn 5 
Co=Key, and Cy=mcy,,* then it is easy to show that 
(i) k=—1 and m==+21 if 

(a) A+C, D+F, K—G, or H+J #0, 
(b) B, LE, or J—H#0, 
(c) A—C=Bm, D—F=Enm, and 2(K+G) =(J—H)m. 





* The results desired here allow ¢,. to be different from zero, and no attempt is made to discuss the 
possibilities in case ¢,—0. 
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(ii) K=—1 and m=0 if 

(a) A, D, H, or G#0, 

(b) B or E#0, 

(c) C=A, F=D, J=H, and G+K=0. 
(iii) K=—1 and m is arbitrary if 

(a) A, D, H, or G+#0, 

(ob). 8=F=—0, C—A, £ =D, J =H, and G+-K—0, 
(iv) k=1 and m=0 if 

(a) A?+B’, D?+E’, or G’+H’+0, 

(b) C=—A, F=—D, K=G, and J=—H. 

The following equations, constructed on the basis of this information, 
answer interesting questions. The equation 

3f[(1+2i)a-+ (3+é)y]—f[(1+3i) 2+ (6+3i) y] —5f (2) +15f(y) =0 
has the general continuous solution f(x”) =¢.(uw?+v’) showing that an equation 
may have a continuous solution, but no non-trivial analytic solution. 

The general solution of 

(1-44) f[3e+ (1—i)y] —3f[2+ 2y] — (6-491) f (x) + (10+ 24) f(y) =0, 
continuous over the finite complex x-plane is ¢.)(u?+2iuv—v’) =cy2, showing 
that an equation may have its analytic solution as the most.general continuous 
solution. 

The elimination outlined in the footnote of $2 shows that every solution 
of an equation of form (1) having no a equal to zero is included in the corre- 
sponding solution of the normal equation whose order is equal to the number 
of distinct ratios 6;/a;, i=1, 2,...., , in the equation of form (1). One 
case remains, namely, that in which each ratio @;/a,;, a; and 6; different from 
zero, is equal to at least one other such ratio and at least one « and one £, of 
different subscripts, are zero. That there are equations of this exceptional 
type which have infinite series solutions is proved by the equation 


f(a+y) —f (ia+iy) —f (—ix) —f(—iy) +f(#) +f(y)=0, i= V—1, 
which is satisfied by a series in positive integral powers of 2* having arbitrary 
coefficients. 


§7. The Converse Theorem. 


It was shown in $6 that any polynomial in x of degree m, satisfies the 
normal equation of order m;. It will now be proved that any polynomial p (2) 
in x, of degree m,, satisfies an equation (1) whose order » is not greater than 
the number of non-vanishing terms of p(x), plus the sum of the degrees of 




















ve 
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such terms, and whose a’s, @’s and y,,, and y,,. may be assigned at will, pro- 
vided a certain determinant A of the a’s and £’s is not zero as a consequence. 


Suppose that 
p(x) =a,e™"+a.0™ + tees +aa™, 


The substitution of p(x) for f(z) in (1) gives an identity from which 
BaMy=—ynur Dar Bly. =0, LBPV=—yeye b=1,2,....,m—1, (12) 
{=1 é=1 i=1 


for the values h=1, 2,...., 7. The total number of independent equations is 
not greater than j+/_,m;. Setting n equal to the number of independent 
equations, we have a system of non-homogeneous linear equations in » un- 
knowns 71, ¥2,---+) Yn,» provided y,,, and y,,, are not both assigned equal to 
zero. A necessary and sufficient condition that the unknown y’s are uniquely 
determined in terms of the a’s, 0’s and two assigned y’s is that the determinant 
A of the coefficients be different from zero. It is obvious that A is a poly- 
nomial in the a’s and @’s. Furthermore, it is at once evident that the term 
formed by the product of the elements in the principal diagonal is unique. 
Hence the polynomial in the «’s and @’s does not vanish identically. There- 
fore the a’s and #’s may be assigned in any way such that the polynomial A 
has a value different from zero. 

If y,4, and y,,,. are both assigned equal to zero, equations (12) form a sys- 
tem of 2 linear homogeneous equations in » unknowns, and the non-vanishing 
of A is a necessary and sufficient condition that each of the unknown y’s is 
zero. In this case the equation is trivially satisfied by p(x). In general, 
then, any polynomial p(x) satisfies an infinity of equations (1) whose orders 
do not exceed the number of non-vanishing terms of p(x) plus the sum of the 
degrees of such terms and whose a’s, B’s, ¥n41’8 ANd Yq49’s may be assigned at 
will, provided only that a polynomial A of the a’s and B’s does not vanish for 
the assigned values. 


$8. Discontinuous Solutions. 


Any solution f(z) of an equation (1) satisfies the normal equation whose 
order n is determined by (1). Suppose that the domain of f(z) is any line in 
the finite complex x-plane, and that f(x”) is continuous in some interval of 
length 6>0 of the line. It is readily seen from the normal equation satisfied 
by f that if the first +1 arguments are so chosen that they represent points 
in the interval while the remaining argument represents a point outside the 
interval, f(a) is determined at the last-named point as the sum of continuous 
functions. Therefore / is continuous at the outside point. In this way f 
may be shown to be continuous at all points in the two intervals of length 
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 6/n which lie at the ends of the given interval. Therefore f is continuous 
in an interval of length [(m+2)d]/n. Any finite interval of length « may 
be reached in this manner by a finite number of extensions of the interval of 
length 6. We may therefore state that if f(a) has a finite point of discon- 
tinuity on any line in the complex x-plane, it has a point of discontinuity in 
every interval of the line, however small. 

Suppose f(z) is continuous in a region of the finite complex w-plane. A 
circle may be inscribed in the region such that f(x) is continuous in the closed 
region of which the circle is the boundary. Suppose the radius of this circle 
is 6 and consider a concentric circle of radius [(n+1)8]/n. By means of the 
normal equation f may be determined at any point in the area between the 
circles as the sum of +1 continuous functions, namely, f at n+1 points in 
the circle of radius 6. Therefore f at the point between the circles is con- 
tinuous. This is true for every point of the area between the circles, and 
hence f is continuous in the circle of radius [(n+1)6]/n. This process may 
obviously be repeated to prove that f is continuous in any finite region of the 
plane. Hence, if f(x) has a point of discontinuity in the finite complex x-plane, 
it has a point of discontinuity in every finite region of the plane. 

G. Hamel (loc. cit.) has exhibited a discontinuous solution* f(x) of the 


Cauchy equation f(a+y) =f (x) +f(y). 


e 


From the treatment in § 2 it is clear that f also satisfies 

f(2u+y)—2f (x+y) +f(y) =0. 
Replacing y by hx+y, multiplying the equation by (—1)’(n—1)!/h!(n—h—1)! 
for successive values h=0,1,....,2—1, and adding the n equations so formed, 
it is easily seen that 


mst. (m1)! 
2 ( 1) 


ki(n+1—k)! 
We therefore have a discontinuous solution of the normal equation for each 
order n. 
§9. Certain Types of Equations having Variable Coefficients. 
The functional equations that have been discussed may be employed to 
solve certain equations of the form, 


2 9.(2, y)f(a2+B:y) +On41(%, Y) (2) +On40(2 Y) f(y) +O,43(2, y) =0, (13) 


where the ¢’s are known functions. <A general statement and a few examples 
suffice to indicate some of the equations that may be solved. Suppose there 
are k transformations 


f(ka+y) =0. (6) 





a=ye'+ey’, y=Aa'tuy’, 





* This solution is obtained on the assumption of the validity of the Zermelo axiom. 









































Witson: On a Certain General Class of Functional Equations. 279 


which may be applied to (13) to obtain new equations such that if each equa- 
tion is multiplied by a non-zero constant, the sum of them is of form (1). A 
solution of (13), if it exists, is included in the corresponding solution of the 
auxiliary equation of form (1) deduced from (13). In order to find a solution 
of (13) it is sufficient to substitute the solution of the auxiliary equation and 
compute the coefficients of the variables. 

Equation (13) includes the non-homogeneous equation in which ¢,(2, y), 
i=1, 2,...., +2, is further restricted to be a constant. In this case it is 
obvious from equation (13) that 9,,3(%,y) is a polynomial in # and y if an 
analytic solution exists, and a polynomial in u, v, s and ¢ if a continuous solu- 
tion exists. Such a non-homogeneous equation is * 


f(x+y)=f(") +f (y) +2ay. (14) 


Transformations which may be used to solve this equation are 
c=2', y=a2'’—y’ and 2=2', y=y'—2’ 
whence, after dropping the primes, 


f (2a—y) —f(x—y) —f(—a#+y) —2f (x) +f(y) =0. 

_ The arguments x—y and —2+y are proportional and the normal equation is 
therefore of order 2. Hence the general solution of (14), analytic over the 
finite complex w-plane is readily seen to be f(z) =a,x+ 2’, where a, is arbitrary. 
The general solution of (14) continuous over the finite complex z-plane is 
f (w%) =ayu+ayv+2’, where a, and a, are arbitrary. 

Suppose all the @’s are constant and ¢,,,—£0. If Ltig,-0, then f(0) is 
finite and uniquely determined, and the transformation f(x) =g(x)+/(0) may 
be employed to obtain an equation (1) of order m in. g(x). Hence f(x) is a 
polynomial (in 2, in u and v, in uw or in v, as the case may be) of degree not 
greater than n. The transformation used in this case has the advantage of 
furnishing an auxiliary equation (1) whose order is not greater than that of 
the original equation (13). 

An equation which illustrates reduction by interchanging arguments of f is 

(cosa) f (x-+y) + (sin?) f (w—y) —f («) —f (y) —2 (cos 2a) ay=0. 
If y is replaced by —y the equation becomes 
(sin’ x) f (x+y) + (cos’ x) f (x—y) —f (x) —f(—y) +2 (cos 2”) ry=0. 
The sum of these equations is of form (1). The solution analytic over the 
finite complex «-plane of the equation having variable coefficients is f(x) =a”. 


It is easily seen that this is the most general continuous solution. 
The equation 


2f(2a+y) + (a+y)f(a—y) +3f(«) —37 (y) =0 


* American Mathematical Monthly, Vol. XXIV (1917), p. 178. 
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may be reduced by replacing x by 2a+y and y by «+ 2y, whence 
2f(Sa+4y) +3 (x+y) f(x—y) +3f(2%+y) —3f(x+2y) =0, 
and subtracting 3 times the original equation from the transformed equation. 
The reduced equation is 
2f(Sa+4y) —3f(2e+y) —3f(%+2y) —9f(x) +9f(y) =0. 
It is readily seen that the equation with variable coefficient has no continuous 
solution. 

If one of the first n+ 2’s of (13) is variable while the remaining ¢’s are 
constant the product of that ¢ and the corresponding f(a+(,y) is a poly- 
nomial. The variable ¢ is therefore a rational function whose denominator is 
a factor of f(a+ 6,y). An equation illustrating this point is 

| S+YT Fey) —1a-+2u) f(a) +2/(y) =0. 
The equation obtained by replacing x by 27+ y and y by w+ 2y is 
2 
18[=*2 | f(x—y) —f (40-+5y) —f (2e-+-y) +2f (2-+2y) =0. 


The equation obtained by subtracting the transformed equation from 9 times 
the original one gives an equation (1) of order 3. The analytic solution f(z) 
of the equation with a variable coefficient is then easily seen to be az’, where a 
is arbitrary. 
$10. Application to Binomial Equations. 
An application of linear functional equations having constant coefficients 
may also be made to certain equations of the form, 


k n+1 
TL [f(ae+By)]"=C Tl (f(ae+Biy) "Lf (y) (15) 


where C is a constant, the real part of each y is positive and no a is zero. 
Let us first consider the solution f(x) of (15) analytic at all points in the 
finite complex plane. Suppose that f(z) has a zero at some point «=a. Let 
y=a in (15). Then k 
TI [f(a@+B,a)]"=0 


for all values of x. Hence there is a finite region in which f(z) has an infinity 
of zeros. But this is impossible since f(x) is analytic throughout the finite 
plane. Therefore, when f(x) is analytic throughout the finite complex plane, 
and not identically zero, it is never zero. The case of a continuous solution 
f(x) of (15) presents more difficulty. It is evident that if one member of (15) 
has either no factor or only one factor involving f, then f(x) is never zero 
unless it is identically so. 
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The function (x) =log f(x), where it is understood that the principal 
determination of the logarithm is employed, is analytic or continuous with f 
when the latter has no zeros. Therefore, in each of the cases considered 
above, @(x) satisfies the equation 


k n+1 
Lye (a,0+B,y) — 2 Ye (a¢+B:Y) —Yn42(y) —K+2sm=0, 


where K is the principal determination of log C and s is any integer. For 
any given s ¢(Z) is a polynomial of degree not greater than n. Furthermore, 
the argumentation of §9 shows that a variation in s affects only the constant 
term of ¢(x). Therefore, 1 each of the cases considered above, f(x) is an 
exponential function of the form 

f(x) =eP**®, (16) 
where P is a polynomial of degree not greater than n, and k(s) is a constant 
depending on s. Thus we see that in general the solutions of (15) are given 
by (16) for the various possible values of k(s). 

The same result may be stated for the continuous solution f(x) of (15), 
when each member involves at least two factors containing f, provided 
{(0)+0. For $(x) as defined, is continuous in some region about the point 
x=0 since f(x) is different from zero in some such region. For a given s, 
therefore, @(z) must be everywhere continuous in the finite complex plane 
because it satisfies a non-exceptional equation of type (1). 

The equation, 

V(a+y)¥(e—y) =[d(2)d(y) 1’, 
mentioned in $1, is included in the last case considered. For suppose there is 
at least one point «=b at which Y(a#) is not zero. Letx=y=b. Then 


(2b) Y(0) =[¥(b) ]*0, 
and ¥(0)#0. It is easily seen that the solutions ~(2) analytic over the 


finite complex plane are 
(x) — et sai —+ ge. 


and the solutions ¥(x) continuous over the finite complex plane are 


W(x) — EC20t? + Cy1UY + Cogd?— srt + Cera. 
where a and the c’s are arbitrary constants. 


§11. Equations Involving More than One Function. 


Consider the equation 


Dyfiaet+ Bey) + ntl n41(%) +Yn+2hns2(y) sn, (17) 
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where no y is zero and the /’s are unknown, continuous, single-valued functions 
to be determined if possible so that (17) shall be identically satisfied by them. 
The functions f; may or may not all be distinct. The method of elimination 
employed in $2 is applicable to (17). If no @ is zero it is evident then that 
f.42 satisfies the normal equation of order n. Therefore f,,. is a polynomial 
of degree not greater than ». Under the assumption that no a; and no £; are 
zero, and no two of the ratios P;/a; are equal, any term may be given the 
argument y by a linear transformation which makes no a and no @ zero. In 
this case, therefore, every function f; of (17) 1s a polynomial of degree not 
greater than n. To find any necessary restrictions on the coefficients of these 
polynomials, it is sufficient to substitute the n-th degree polynomials having 
general coefficients in (17), and to equate to zero the resulting coefficients 
of the variables. It is obvious that equation (1) is a special case of (17). 

If some of the ratios 8;/a; are equal it may be assumed without loss of 
generality that equation (17) is so arranged that functions having arguments 
of a common ratio are placed consecutively. If all the functions having argu- 
ments of a common ratio have subscripts 1, such that geich, then we may 


write 


F,(a,0+ By) = 2h (a~+ Py). 


Equation (17) may now be written 

LF, (a,0+ Bry) +Fagi(%) +Fis2(y) =9, (18) 
where no a and no @ are zero, and no two ratios 6,/a, are equal. We denote 
by g+2 the number of terms in the first member of (18). Each F is a poly- 
nomial of degree not greater than g. Each F therefore determines a non- 
homogeneous equation in certain f’s having arguments differing by constant 
factors. If the f’s of any F are identical, that is, if the f’s of F are the same 
function, F determines a non-homogeneous mixed g-difference equation satis- 
fied by f. The equations of type (1), which have no a@ equal to zero, but have 
some ratios B;/a; equal, or some $’s zero, are special cases of (18) given by 
Fis2(Y) =Yns2f(y). The equations of the exceptional case noted in the last 
paragraph of §6 are equations of form (18) which have no F a constant 
multiple of a single f. In this connection it is interesting to note that the 
function f of the example in the paragraph cited satisfies the equations 


f(«)—f(—tx) =f(«) —f (ix) =0, 
f(—tw) =f (tx) =f(—a#) =f (2). 


whence 



































Contributions to the Study of Oscillation Properties of the 
Solutions of Linear Differential Equations 


of the Second Order." 
By R. G. D. Ricuarpson. 





Introduction. 


The study of boundary problems for linear differential equations of the 
second order dates back to the time of Euler and D’Alembert, with whom it 
arose in connection with problems of mathematical physics. Beginning with 
the fundamental paper of Sturm in 1836, there have been extensive investiga- 
tions ¢ in this field in recent years, notably by Klein, Bécher, Stekeloff, Kneser, 
Hilbert and Birkhoff. Since the differential equation of the second order is 
of such fundamental importance in so many fields, and since similar general 
problems for equations of higher order can not be handled by processes so far 
devised, the invention of new methods and further investigation of the nature 
of solutions find ready justification. 

The equation to be studied will be taken in the form 


d dy (x 

f (w(a) Y)) 4 G(x, a)y(x)—=(py.).+G(m, X)y=0, OSeS1, (1) 
where G(x, A) is a function depending on a parameter a. The solution y(z) 
of this self-adjoint equation shall be subject to the self-adjoint boundary con- 


i my (0) +anye(0) aay (1) -+ax¥%e(1) =0, 
B1y (0) + Bry. (0) +Bsy (1) +B.y.(1) =0, (2) 
p(1) (a82.—%81) =p (0) (a8,—a,83), 
where the two sets of real coefficients a,, a,, a,a, and 6,, B,, Bs, By are 
linearly independent. The most important special cases of these boundary 
conditions are given by 
y(0)=y(1)=0; y(0)=y,(1)=0; y,(0)=y(1)=0; y,(0)=y,(1)=0. (2’) 








* Read before the American Mathematical Society, September 4, 1917. 

} For existing methods and literature of the subject see Bécher, Encyklopidie Mathematischen 
Wissenschaften, II A7a, Proceedings International Congress of Mathematicians, Vol. I (1912), p. 163, and 
* Lecons sur les Methodes de Sturm” (1917); Lichtenstein, Rendiconti del Circolo Matematico di Palermo, 
Vol. XXXVIII, p. 113. 
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There are many interesting questions in regard to this linear problem. Do 
there exist parameter values 4 such that there are solutions of (1) satisfying 
relations (2)? If so, how many are there, and how are they distributed? What 
is the nature of the corresponding solutions? How do the solutions vary with 
change of the coefficients of the equation and of the boundary conditions? 
When does the totality of solutions form a fundamental set in terms of which 
functions may be expanded? 

-Among the methods which have been used in studying the problem are: 
(1) Differential equations including the use of comparison, approximation and 
asymptotic expressions; (2) the minimum principle in the calculus of varia- 
tions; (3) integral equations; (4) the theory of linear algebraic equations 
in an infinite number of variables;* (5) a limiting process with linear 
algebraic approximating difference equations (cf. §1, IV). The methods 
developed in this memoir would fall under (1) and (2) and may be character- 
ized as a free use of the differentiation of fundamental formulae with regard 
to the parameters involved in the equations and in the boundary conditions, 
together with the proof that under the conditions imposed certain integrals 
are positive. 

Exact oscillation theorems for solutions under the boundary conditions 
(2) have been developed by Birkhoff t+ for the case of the special equation 


ee ee ee a a 


OA A=—0 A=+0 


Another important special case of equation (1) in which G(z, 4) contains the 
parameter linearly 


(py:)2+ (Qq+ak)y=0, (3) 


has been studied very extensively. When k(z)>0 this equation may be 
reduced to a form included in that investigated by Birkhoff. 
The definite case of (3), viz., when one of the integrals 


1 1 
f ky'dx, f (vyi—ay*)de 
0 0 


has one sign for all functions y(z) considered, has been discussed in many 
phases by mathematicians since the time of Sturm. By means of his theory 
of integral equations Hilbert { established the existence of characteristic 








* Lichtenstein, loc. cit. 


+ Transactions of the American Mathematical Society, Vol. X (1909), p. 259. It should be remarked 
that both in the second and third lines from the end of the statement of the principal theorem (p. 269) 


instead of p-++ 1 we should read p—1. 
¢ “Grundziige einer allgemeinen Theorie der linearen Integralgleichungen” (Teubner, 1912). 
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parameter values and characteristic solutions of equation (3) for this definite 
case. As will be shown in $5 the boundary conditions which he used are nor- 
mal forms of (2). Oscillation theorems for equation (3), under some of the 
simple boundary conditions (2’), have been established by various means; 
among others by setting up the corresponding calculus of variations problem, 
and interpreting the Jacobi criterion.* 

When k has both signs and q is positive in at least a part of the interval 
and sufficiently large there, it is not necessary that either of the integrals in 
question be definite. This case, which we shall call the non-definite, was first 
discussed incidentally by the author in a papert in which he was treating the 
problem of oscillation theorems for two equations with two parameters. He 
showed that when the boundary conditions are y(0) =y(1) =0, there exists an 
integer n, such that for n<vn, there are no real solutions which have n zeros, 
while for n>, there are at least two.t At that time all the principal results 
of $$ 2-4 were obtained, but were not published. 

For the special equation (3) and k>0 the theorems proyed by Birkhoff 
were rediscovered by Haupt§ in his dissertation. The oscillation theorem 
stated in this dissertation for the case that k changes sign is corrected in a 
later article,|| and by using the methods which I had developed, various oscil- 
lation theorems for the general equation (1) are derived. It is also shown by 
means of expansion theorems that if the non-definite equation (3) be taken in 
a certain normal form there exists an integer n, such that for n>n, there are 
precisely two solutions with n zeros and satisfying the boundary conditions (2). 

The object of the present memoir is to investigate the conditions to be 
imposed on G(x, A) regarded as a function of A, so that definite oscillation 
theorems for solution of (1) may be determined. In § 2 an attempt is made 
to bring as close together as possible necessary conditions and sufficient condi- 
tions for a limited or unlimited number of oscillations. A criterion for the 
behavior of the zeros with change of the parameter is obtained in §3. This 
permits the development of very general theorems for the unique existence of 
solutions vanishing at the end points and possessing a prescribed number of 
zeros, and also for the existence of two, and only two solutions of this nature. 





* Mathematische Annalen, Vol. LXVIII (1910), p. 279. 

+ Transactions of the American Mathematical Society, Vol. XIII (1912), p. 22. 

+ That there were exactly two when n>n, was stated by the author ina paper in the Mathematische 
Annalen, Vol. LXXIII, p. 289, in which he was discussing oscillation theorems for three linear equations 
with three parameters. The subsequent results of the memoir were not affected by this error which was 
corrected by a note in Vol. LX XIV (1913), p. 312, of the same journal. 

§ “ Untersuchungen iiber Oszillationstheoreme ” (Teubner, 1911). 
|| Haupt, Mathematische Annalen, Vol. LXXVI, p. 67. 
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These theorems contain as special cases all known results in this field, and 
some new special cases are set forth in detail. 

The non-definite case of (3) and (1) is discussed in §4. The question of 
whether there may be for a given oscillation number more than two parameter 
values for which there are solutions of (3) (2) is settled by giving an example 
in which for any n in an interval n,, n, there are four values of A corre- 
sponding. That there exists an integer nm, such that for n>n,, there are 
exactly two solutions is proved by a method entirely different from that of 
Haupt, and in some particulars it would seem that the resulting theorem is 
less satisfactory, in others more satisfactory than his. The theory is also 
extended to cover some corresponding cases of (1). Concerning the complex 
solutions which correspond to values of n<n, some theorems are derived. 

In the latter half of the memoir a method is developed for obtaining the 
facts in regard to the solutions of the general equation (1) under the geueral 
boundary conditions (2). With this end in view $5 is devoted to a reduction 
of the boundary conditions to normal forms by means of the usual transfor- 
mation of the dependent variable y, which leaves the number of zeros 
unchanged. Each of these three normal forms, 


I. cy(0)+y,(0)=0, ty(1)+y,(1)=0; 
II. y (0) =hy(1), hp(0)y,(0) =p(1)y,(1); 
III. y(0)=lp(1)y,(0), lp (0) y,(0) =y(1) 


contains one or two parameters o,t, h,/. When special values 0 or o are 
assigned to these parameters, the forms reduce to the simple cases (2’) for 
which the facts are readily obtainable from the developments of the earlier 
sections. The first form is of essentially different character from the others. 
In the latter the parameter h or / is a double-valued function of A, which 
in general is real in sub-intervals only, while in the former each of the 
parameters c, t is a single-valued function of the other and of A, and is real 
throughout. By letting 4 vary, and calculating the rates of change of these 
parameters and of G(z, 2), theorems of oscillation are derived (§{$6-8) for 
boundary conditions in each of the normal forms. While detailed results are 
not given in all cases, these are immediate developments of the fundamental 
facts ascertained. 
Any linear differential equation of the second order 


PYe2tVy2tOy=0 or py. +by.+ (A%+6)y=0, o>0, 
may be thrown into the corresponding self-adjoint form (1) or (3) on multi- 


"Yaz ; eels 
oo" ; the corresponding self-adjoint boundary 


plying by the function p(z) = 
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conditions (2) are changed in form only by the substitution p(0)=1, 
6% dz 
p(l)=err? 

The oscillation theorems remain unchanged if the variables are subjected 

to the usual transformations of the dependent and independent variables 


y=ny, w= (£) (n(2) +0, “ #0). The resulting equation is of the form 


PVc2t+ Vet 07=0 or FHiz+V9z+ (A, +6,)9=0, 

the new boundary conditions for the new interval Z, 2, being of the form 

a1Y (Zo) + ais (Zo) + asY (%1) +49; (1) =0, 

BxY (Xo) +BY (Zo) +BsY (71) + Bz (ts) =9, 

oF 2 - ei == = Ss 
Je 6 [a,0,.—428,] =4,0,—a,Bs, 

as may be shown by computation.* The invariantive property of self-adjoint- 
ness gives to the results obtained for (1) or (3) a very general character. 


$1. Some Properties of Solutions of the Differential Equation. 
In his original memoir Sturm studied the differential equation in the form 
d = ai =~ 


where K>0 and G both depend on a parameter 2. But by a change of 
variables this may be reduced to the form 


(pyz)2+G(a, A) y=, (4) 
where p(x) is a positive function independent of a. We lose nothing in 


generality by considering this latter equation. Regarded as functions of x 


the coefficients p, G will be postulated as continuous together with as many 
derivatives as is desired, while G will be considered as analytic with respect 
to A. The usual modifications of the results derived can be written down im- 
mediately if less stringent hypotheses are imposed. The trivial solution y=0 
will be excluded from the discussion. Some theorems concerning solutions of 
(4) will now be reviewed. 

I. If one zero of a solution of (4) is held fixed, all others are moved 
nearer to it by a decrease of p or an increase of G. If, for example, p(x) is 
less than a constant P, and G(a, A) is greater than a constant y>0, the zeros 
of (4) are closer together than those of the equation 


Yort gy=0, I= z ?. 





*For the special case of a transformation of the dependent variable only, this computation is 
given in §5. 
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which has a solution y=sin Vg(x+c) with zeros at intervals of Vi . By 


fixing p, and taking G large enough in any interval of x, the zeros of (4) may 
then be made as close together as desired. 
II. The special equation where G(x, A) contains the parameter linearly, 


(pyz)2+ (q+ak)y=0, (5) 
has been much discussed. The boundary conditions 
y(0)=y(1) =0 (6) 


are of the greatest interest and three cases may be distinguished. 

(A) Orthogonal Case, when k(x)20. There is then an infinite number 
of parameter values 4,,(A,<A,<....) with a limiting point at positive infinity 
only, for each of which a solution Y,, satisfying (6) exists. The number of 
zeros of the soluiion Y,, (including those at z=0 and x=1) is m+1. 

(B) Polar Case, when k(x) takes on both signs and the integral 


1 
D(y)= J (pys—ay*) de (7) 
is positive-definite,* that is, for the given boundary condition (6) D(y) can 
not take on negative values. There are two sets each of an infinite number of 
parameter values 0<2,<”%)...., 0O2”2_,24_,...., with limiting points at 
positive and negative infinity respectively, corresponding to which solutions 


Y,,, Y_n exist. For both Y,, and Y_,, the number of zeros is m+1. 
(C) Non-definite Case, when both of the integrals 


1 5 
f hyd, § (pyt—ay*)da 
0 0 
may take on negative values. This will be discussed in § 4. 
III. For the equation (5) and the special boundary conditions (6) cer- 
tain minimum properties may be stated. In the orthogonal case the minimum 


of the integral D(y) for those values of y which satisfy (6), and the normal- 
izing and orthogonalizing conditions 


1 1 
f kyde=1, f k¥.ydx=0, i=1, 2,....,m—1, (8) 
0 0 


is A,,, and is furnished by the normalized solution Y,, of (5). In the polar 
case the minimum of D(y) for those values which satisfy conditions (6) and 
(8) is A,, and is furnished by Y,,; the minimum subject to the conditions (6) 


and 1 1 
f ky*dx=—1, f kY_.yda=0,. +=1, 2,....,m—l, 
0 0 





* Bécher has pointed out (Proceedings International Congress of Mathematics, loc. cit., p.173) that 
the special case of the polar problem where g <0 can be treated by the method of Sturm. This remark, 


however, does not apply to the polar case in its most general form. 
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is —A_, and is furnished by Y_,,. In the non-definite case, D(y) can be nega- 
tive and the minimum (even in the simplest problem (m=1)) may not exist. 
However, in view of the developments of § 4, it would seem probable that by a 
modification of the discussion, the solution Y,,, for m large enough, may be 
regarded as furnishing a minimum of a calculus of variations problem. 


IV. While the notion of regarding a differential equation directly as the 
limit of a set of difference equations has been used heuristically since the time 
of Euler, it was, so far as the author is aware, first made definite and rigorous 
in a recent paper.* In that paper the problem actually discussed is that of 
existence theorems for partial differential equations with given boundary con- 
ditions. But the same method applies, for example, to equation (1), and the 
discussion is essentially simpler for the equation in one dimension than for 
that in two or more. For the sake of simplicity let us confine ourselves to the 
case (5) (6) and consider the unit interval to be divided into m equal parts, 


the values of y, p, g, k at the point —(i=0, 1,....,m) to be denoted by y;, p;, 


di, k;, and difference equations 
m*[ Di+1 (Yis1—Y) —p:(y¥i—Yi-a) | tay:tAaky,=0 (1=1,2,....,m—1), (9) 


to be set up. In order that there be solutions of these equations ~ must be 
one of the m—1 roots of the determinant formed from the coefficients. With 
increase of m the number of points at which y is defined increases, but we can 
pick out corresponding parameter values and solutions of the various sets of 
difference equations, and if proper continuity conditions are imposed on the 
coefficients of the differential equations, it may be shown that the correspond- 
ing sets of parameter values approach as a limit a parameter value of (5), and 
corresponding solutions approach a solution of (5). In this way the infinite 
set of solutions of the differential equation is obtained. If & has both signs 
and q is positive and sufficiently large, at least in some part of the interval, 
some of the A’s and the corresponding y’s will be complex.t But in all those 
cases of equation (5) heretofore treated (the orthogonal and polar), the 
method of passing to the limit in (9) suffices and gives a simple proof of the 
fundamental facts. This method can be extended to a treatment of existence 
theorems for solutions of the equation (4) with more general boundary con- 
ditions. 





* Transactions of the American Mathematical Society, Vol. XVIII (1917), p. 489. 
+ One method of proving this would be by noting that the parameter values and solutions are 
approximations to those of the differential equation which are shown in §4 to be complex. 
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§2. Sufficient Conditions for the Existence of Solutions with an Unlimited 
Number of Zeros. 
In the consideration of the equation 
(PYz)2+G (x, A)y=0 (10) 
let us impose the restriction that for finite values of 4 the function G(2, A) is 
limited, and in particular that G(2,0) is limited. The interesting cases will 
be covered by one of two hypotheses, which will be justified by the later devel- 
opments of this section. 
Hypotuesis A. For at least one point of the interval, the upper limit of 


G(x,%) becomes infinite with a (lim G(a,2)=+0), and in such a manner 
A=+0 


that the number of zeros of the solutions increases without limit with a. 
The problem treated by Birkhoff* where a >0 and G(a, +o0)=+o is 


a special case, and the orthogonal problem (§1, II) is still more special. 
Hyporuesis B. For at least one point of the interval lim G(%,a)=+o0; 
A=+0 


for at least one other, lim G(a,2)=-+0; and in both cases G(x, 2) increases 


A=—-2 


in such a manner that the number of zeros of the solutions increases without 
limit with 2. 

The polar case (§1, II) is included in this hypothesis. 

THEOREM I. In order that there be a set of parameter values A such that 
the number of oscillations of the corresponding solutions of (10) be unlimited, 
it is necessary that in the neighborhood of at least one point lim G(a,A) =+ 

jects aes : A=+o0 
or im G(x, A)=+o. 

For, as the number of zeros increases, the length of the smallest interval 
decreases without limit. The zeros of the equation 

Py,,+G(a,a)y=0, P=maximum p(z) (11) 
are farther apart than those of (10) ($1, 1). To establish the theorem for 
(10) it is then only necessary to prove it for (11). Let us denote by a, @ that 
pair of consecutive zeros of (11) whose distance is a minimum. In such an 
interval a, 8 [y(a) =y(8) =0], y may be taken positive, and since the equa- 
tion is homogeneous, its solution for all values of 4 may be multiplied by a 
constant so that the maximum is 1. The maximum of y, must be at least as 


great as 7 which is the slope of the line joining (a,0), (8,1). To investi- 





* Loc. cit. 
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Yuu 





gate the maximum of — re we note that y, is zero at some point of the inter- 
val and at least as great as I at another. Hence —y,, must be at least as 
1 P y 1 ‘ 
reat as —.——, at some point, and further max (— ua) > ———.,. Since 
. (B—a)? — y / = (B—a)? 


for at least one sub-interval the length approaches zero, and since 


P Ee 7 

G(a,a)=— - . 
it follows that for the interval 0, 1 lim max G(x, 2)=-+. Since by hypothesis 
G can become infinite only for A=+o the theorem may be readily deduced 

from these results by the usual processes of reasoning. 
To show that the necessary condition of Theorem I is not sufficient, let us 
consider the following example: In the interval 0, 4—e we set up the function 
9 

y= = sin, and in the interval $+<¢,1 the function y= = sin (1—2z); 
both of these ares satisfy the equation y,,+ = y=0. If continued, they 


would meet at the point x=4 with an angle arc tan 2\/2. If in the interval 
4—e,4+e an analytic curve is introduced which is tangent to these two arcs 


(on one of which y,> and on the other y,< ), it must have at some 


1 1 
v2 V2 


point —y,,> and if ¢ is taken small enough, max (— vss) is great at 


zs 
V2e’ 
pleasure. Hence denoting by y,,+G(z,2)y=0 the differential equation which 
has for solution the function defined for the interval 0, 1 by this method, and. 
setting a=, we know that lim max G(z,4) in the neighborhood of x=43 


e=0 : 
increases without limit. On the other hand none of the suite of functions y 
has any zeros in the interval. 
Having shown by this example that the necessary condition of Theorem Lis 
not sufficient, let us now deduce a sufficient criterion that when lim G(2,A)=0, 


the number of zeros of solutions of (10) be unlimited. A similar discussion 
can be given for the case lim G(2,4)=o. Denoting as before by P the maxi- 


A=—a 
mum of p(x), and by M an arbitrarily large constant, it is possible to find a a 
and an 2-interval such that G(z,4)>MP in this interval, whose length will be 
denoted by e,. The number of zeros of the solution sin \/M(xz+c) of the 
equation y.,+My=0 in an interval of length ey is not less than the integral 
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part of vee It follows immediately from §1, I that the number of zeros of 


solutions of (10) is not less than that of sin VM(a2-+c). Hence the 


Tuerorem II. In order that there be N zeros in a solution of (10) it is 


‘ M , 
sufficient that one can find an M such that yom >N, where ey is the length of 


a sub-interval where G(2, 2) >MP. 
Corotuary I. In order that the number of zeros be unlimited it is 
sufficient that lim ey VM=o. 


A= 
Corotuary II. If throughout any sub-interval of fixed length the value 
of G(x, A) increases without limit, the number of zeros increases indefinitely. 


However, the sum of the sub-intervals in which G>M may remain above 
a positive constant without compelling the number of zeros to increase with M. 
For example, one can set up a function which has no zeros except at x~=0 and 
x=1 and which, except in the neighborhood of these points, oscillates between 
y= and y=+3, the number of oscillations increasing indefinitely with A. More- 
over, by taking the oscillations frequent enough one can have G(x, 2) >MP in 
portions which total at least one-quarter (or any other proper fraction) of the 
interval, the great curvature downward in these sub-intervals being counter- 
balanced by curvature upward in those where G(a, A) has the opposite sign. 

On the other hand, we have seen that the number of zeros of the solution 
sin \/M(x+c) of the equation y,,+My=0 is not less than the integral part of 


= where e, is the length of the interval]. If this interval is divided into 


two parts ¢,, e, for which there are solutions sin \/M(x+c,), sin \/M(2%+<c,) 
respectively, the number of zeros can not be reduced by more than one. For, 


the sum of the integral parts of = and _ eannot differ by more than 





one from the integral part of ete) . It follows in the same way that if 


the interval is divided into n, parts the number of zeros of the solutions can 
not be decreased by more than v,—1. But in the intervals where G(z,a) >MP 
the zeros of the solutions of (10) must be at least as many as the minimum 
number for y,,+My=0. 

Tueorem III. In order that the solution of (10) have N zeros it is 


sufficient that one can find an M such that vary —n,+1>N, where ny denotes 


the sum of the lengths of the n, intervals in which oe >M. 
max p(2%) 
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§3. Behavior of the Zeros with Monotone Change of i. 


Let us fix the zero at the left-hand end of the interval 0, 1, denote such a 
solution by Y [Y(0)=0], and investigate what happens to the other zeros as 
4 increases. It is convenient to think of G(a,’2) as being defined for values 
x>1. Since the coefficients of (10) are continuous, the zeros move continu- 
ously, and since we can not have both y(z)=0 and y’(x)=0 without having 
y=0, the zeros can not coalesce and then disappear. By differentiating (10) 


with regard to A we have 


O(DYs) 2 
a. +a5t 4 Fy sai, 


and on multiplication of this equation by —y and of (10) by > addition and 


(12) 


integration from 2, to x, we get the fundamental formula * 
oy | | Sue" . aes 
Ba lt ih aa 9 (13) 
. : oY (0) iia 
For the particular solution Y (xz) we have Y(0)=0, a and if a is another 


zero of Y the formula becomes 


oY 204 o» ; 
p(a)¥,(a)5-(a)=f 5 Fide. (13’) . 


Y(a) 


If the integral on the right is positive, the signs of Y,(a) and g aa, are the 


same. Hence when Y,(a) is negative, Y(a) is decreasing with increase of 4, 
and this zero of Y is moving to the left; when Y,(a) is positive, Y(a) is 
increasing and the zero is again moving to the left. In the same way it may 
it may be shown that all zeros for which the integral on the right of (13’) is 
negative, move to the right. When the integral is zero, further investigation 
is needed. We have now proved 





*If p is also a function of \ this formula becomes 


y OY —[ 2 oe . tz Op y? 
PYz S¢ | ae |. Sxueu | =[-s y*da— a & an’ ‘2a. 


In terms of this hypothesis results wp to all the Mescstine theorems of this paper may be written 
down. Since, however, by a transformation of variables p can be made independent of A we shall, for 
the sake of simplicity, confine ourselves to the problem proposed. 

Oye(0) | Oy) _ 


If at x=—0 the condition y, (0) + cy (0)—0 is imposed where oc is a constant, then ar +o Da 0, 


and this formula (13) can be written 
Oy >) x 
= Oy _,, OY as . 
vP = (8) =| 1.3% y al =\f\a da. 


Hence, in order that the roots of a solution y and the roots of or alternate it is sufficient that 


zs oe y°dew have one sign. 
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TuroremM IV. If Y(0)=0 and tf a is a parameter value for which there 
are subsequent zeros of Y at a,,....,a, then with increase of A a zero, a,;, 
moves to the left or right according as 


f° SE nw, ayae 


is positive or negative. 
The argument used in deducing Theorem IV is still valid if in place of a 
zero of y we introduce a zero of y,. 


Tueorem IV A. If y(0)=0 or y,(0) =0 and if A is a parameter value 
such that y(a)=0 or y,(a) =0, then with increase of A the zero a moves to the 


left or right according as ; 
°0G(%,A) , 
J a aaa (x, A) dx 
is positive or negative. 
Corotuary. In order that 2 be a multiple characteristic number and 
Y(x,2) a multiple solution of the problem (10) (6) it is necessary* that 


ee ¥*(x,2)de—0. 


G can be so chosen as a function of A that with increasing 4 the integral 


1 0G 
OA. 


wiaaiee of A may be small at pleasure, we see that a zero may pass back and 
forth through x=1 an infinite number of times. There may then be an infinite 
number of A’s for which exist solutions of (10) (6) with a fixed number of 
ZeYros. | 

As A runs through its interval there will be at least one value for which 
the number of zeros of Y(#,4) isa minimum. For both the orthogonal and 
polar cases this minimum number of oscillations is zero, but in general this 
will not be the case. If for A=A”’ there is a minimum number 2, of oscilla- 
tions, then from the principle of continuity we can argue that for n>n, there 
is under Hypothesis A (§ 2) at least one value of 4 for which exists a solution 
of (10) vanishing at <=0 and w=1, and oscillating » times (having n zeros) 
in the interval, while under Hypothesis B there are at least two. 

One fundamental problem for investigation is the determination of 
sufficient conditions that in a given interval of 4 (which may include infinity) 
there is only one solution of (10) (6) oscillating a given number of times; in 


x “dar i is alternately positive and negative, and since the corresponding 





* This condition may be shown to be sufficient also. 
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other words, sufficient conditions that the zeros of Y (a) pass through the point 
x=1 in one direction only as 4 increases (or decreases). Another problem is 
that of determining when there exist precisely two solutions oscillating a given 
number of times. 

Tueorem V. I[f for values of Ain an interval (which may include infinity) 
there are solutions of (10) (6) for which the minimum and maximum number 
of oscillations are n, and n, (nz, may be infinite), respectively, and if the value 


of ns is equal to or greater than o(A)G(a, A), where > is a positive function 
of A, then there is one and only one solution oscillating n times (nySn<nz) m 
the interval 0<a”<1. 

For, we have on multiplying (10) by Y, and integrating under the boun- 


dary conditions (6), 
1 1 
f pY¥%dx= f GY?da. (14) 
0 0 


From the hypothesis we may then write 
1 0G wo " re we 
J 55 Vidw2oa) f GY'da=9(4) f pY¥%da>0 


and since this holds for all A, Theorem IV gives the desired result. 

The theorem just proved fits the case where with unlimited increase of a 
the maximum of G increases without limit. For the case where 4 approaches 
—oo, we observe from Theorem IV that a sufficient condition for a similar 


theorem is that if = Y*daz be negative. To insure this, it is sufficient to take 
- 9 


ag 
OA 


<—G, where (A) is a positive function. Then 


ag - 
5 Vida <—4(2) £ GY*da=—W(A) i} pY2%dx <0. 


Combining this result with Theorem V we may state the following: 


TurorEeM VI. [f for two values 4’, A” (A’ =A") of A there are n, zeros of 


the solution Y within the interval 0,1, and if for A>2’, SG > o,(a) Ga, Xr); 


and forra<a", La. < —,(A)G(a, 2) where $,, o are positive functions defined 


in the intervals 24’, +0 and —ow, A”, respectively, then under Hypothesis B 
(§2) for n2n, there are two and only two values of A for which there is a 
solution Y(x,A) of (10) (6) with n zeros. 
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In the polar problem we have the special case of the hypothesis of this 
0G 1 


theorem where 4’=A”=0, n,=0, G=q-+Ak, >a, =, )= = >0, —9?,= +<0. 
A more general special case of Theorem VI is discussed in § 4. 
Since formula (14) is valid also when 
y(0)=y,(1) =0, or y,(0)=y(1)=0, or y,(0) =y,(1) =9, 
theorems analogous to VI can be written down for any of these boundary 
conditions. 
§4. The Non-Defimte Case. 


With regard to applications, the differential equation 


(DYz)2+ (a+ak)y=0, (15) 
with boundary conditions 
y (0) =y(1) =0, *(16) 
is the most important special type of (10). The problem where one of the 
integrals 


1 ; 1 
D(y, 9,1) = 7 (pyz—qy’) da, J ky*dax 


is definite, has been studied in great detail. Let us first investigate some 
properties of the remaining non-definite case where both these integrals may 
take on negative values and later deduce some analogous results for the more 
general equation (10). In $4 6-8 these results will be extended to cover the 
more general boundary condition (2). 

For qg<0 the integral D(y, 0,1) can not be negative. If qg is positive in 
at least a part of the interval 0, 1, the minimum of the integral’ Sf pyzdxz under 
the conditions (16) and f*gy’dz=1 is (§1, III) the smallest parameter value 





A of the equation (py,),+Agy=0, and when —# pa™ this in turn is ($1, I) 


at least as great as the smallest of the equation y,,+Am’y=0 under the same 
boundary conditions (16). Since for the latter equation 4,=1 and Y,=c sin aa, 
it follows that 4 is greater than unity. Hence f' py2dv=A is greater than 
Sf. ay’dx=1 and D(y, 0, 1) is positive. 

q 
min p 
part of the interval, and this assumption we shall now make. But we note 


here that in a stb-interval a,, a, such that y(a,) =y(a,) =0, any value of q for 
2 


>z’ in at least a 





In the non-definite case one must therefore have 


7 
(@_—a,)*’ 





which —A : < makes the integral {7 (py:—qy’)d« definite, as is 
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readily shown by a comparison argument like that above. Hence, for given p 
and q, if throughout the interval 0,1 the zeros can be taken close enough 
together (as e.g. may certainly be done where A has the same sign as k(x) and 
is taken large enough in absolute value, cf. §1, I), the integral for each sub- 
interval is definite, and hence D(y, 0,1) is positive. This is also true under 
the more general boundary conditions (2). 

On the other hand from the formula 


D(y, a1, %)= J (pyi—ay')da=Af “ky*da, (17) 
which is obtained by multiplying (15) by y and integrating under the condi- 
tions y(a,)=y(a,)=0, we see that if A is positive and k negative in the 
interval, then D(y, a,, a) is negative; it is, however, still possible that for 
the larger interval 0, 1 the integral D(y, 0,1) be positive. Were this the case 
(as will later be established provided A is large enough) we could deduce from 
Theorem IV and the formula 


1 1 
D(y, 0,1) =f (pyt—ay*)da=af ky'dx, y(0)=y(1)=0, (18) 
0 0 


certain facts in regard to the movement of the subsequent zeros of particular 
solutions Y of (15) which vanish at x=0. Before proceeding to a detailed 
discussion of this matter let us compare some possibilities for the orthogonal, 
polar and non-definite cases. 

In the orthogonal case let us consider all 4 greater than some finite num- 
ber chosen less than the smallest characteristic (which may be positive or 
negative). As A increases, the subsequent zeros of any solution Y of (15) 
[Y (0) =0] move to the left, new ones being added to the interval 0,1. For 
the smallest value of 4 there are no zeros within the interval and when A 
passes through a parameter value 4, the n-th zero enters. In the polar case 
there are two ranges of values of 2; for the range which extends from zero to 
+o, a result precisely like that of the orthogonal case may be stated; for the 
range which extends from —o to zero, a decrease of 4 causes subsequent zeros 
to move to the left, there being no zero of Y(x,0) present in the interval while 
the n-th enters for A=A_,. In the non-definite case we know that the results 
must be quite different and we may expect that the march of the zeros will not 
be monotone with 4. In fact, there may be a range of values of 4 (L,5a<L,) 
such that as 2 increases the number of zeros first decreases, then increases, 
then decreases and finally increases, the minimum number being a positive 





integer.” 


* A reference to the proof of this last fact is given in the Introduction. 


38 
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As an example let us consider the equation 
Yoot [(1002)*+Ax)y=0, y(0)=y(1)=0, 

where x is equal —1 in the sub-interval 0, $ and equal 1 in the remainder of 
the interval. Such a function x(x) may be approximated by an analytic func- 
tion k(x) for which the corresponding equation has solutions with similar 
properties. When ~2=0 the solution y=sin 100zz has 100 zeros in the interval. 
When A=(1002)?, we have the equation y,,=0 in 0, 4 and y,,+2(1007)’y=0 
elsewhere, the solutions being respectively y=ca (c=const.), y=sin 100\/27a. 
There are are then approximately 50\/2 zeros. For A=—(1002)* there are 
evidently the same number of zeros. More generally, when 4 is in the interval 
—(1002)*, (1007) the number of zeros is approximately 


= [ V (10072)?+ a+ V (1002)?—a], 








a function which has its maximum for A=0 and decreases with increase of ||. 
On the other hand when |A| is greater than (1007)* the number of zeros is 





approximately = V 00%)" + |A|, which increases with |A|. The minimum 


number of zeros is then approximately 50\/2 and occurs for A= (1002)? and 
—(100z)*.* For values of n between 50\/2 and 100 there are four parameters 
corresponding, to which exist solutions oscillating n times. 
Returning now to the discussion of the sign of D(y, 0, 1) we can prove 
Lemma l. If in an interval a, b, the function k is positive, except perhaps 
at the end points, then 2 can be taken so large that for solutions of (15), 


b 
D(y, a, b)= { (pyi—ay*)da>0. 


For, by taking A large enough, q+Ak can be made as large as desired, 
except perhaps at the end points. Hence the zeros are as thickly strewn 
throughout the interval a, b as desired ($1,1). We note that in any sub- 
interval 8,, 8. of a, b, in which there is a zero of y and for which the maximum 


of jy| is M, the minimum value of es y2dx is taken on when y is a straight lize 
f2 


joining (@,,0) and (6,, +M). This minimum value is a a Since on the 
as cs | 
other hand Spayda< (C.—P1)M? max q, the value of D(y, @;, 0.) is greater 


than | es — (P,—,) max a|: hence when (8,—,)* is made less than 
2 1 


* By reference to Theorem IV it follows from this discussion that in the two intervals which are 
approximately —(1007)*,0; (1007)*, 4-0 the value of Siayde must be positive or Siyda<fiyda, while 
in the: intervals which are approximately —o,—(1007)?; 0, (1007)°, Lb y'da>f y’da. 
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= this integral is definite. It follows that since the first and last zeros of 
y(x) are as close to the ends of the interval a, b as is desired, those portions 
of the integral arising from these two end sub-intervals can be made positive. 
For any interval a,, a, of oscillation [y(a,)=y(a,)=0] it follows at once 
from (17) that when k and 4 are positive D(y, a,,a,) >0. Combining these 
results we have the lemma. 

It follows in the same way that if k<0 and k$0 in an interval a, b, then 
2 can be taken so large negatively that for solutions of (15) the integral 
D(y, a, b) is positive. 

Coming now to an interval c, d in which k is negative, except perhaps at 
c and d, it is possible in an interval c+y, d—y (y>0 arbitrarily small) to 
take M’ arbitrarily large and then choose A so large that —(q+Ak) >M’ or 


| (pye)2|=|—(at+ak)y|>M’ly|, [(py.)2>My if y>0]. = (19) 
The are y(x) may be shown to be sharply concave away from the z-axis.” It 
is possible to prove the following: 








Lemna II. If in an interval c, d the function k 1s negative, except per- 
haps at the end points, then 2 can be taken so large that for a solution of (15) 


ad 
D(y, ¢, d) = { (py:—aqy’)da>0. 


This will be proved by showing that on taking 4 great enough it is 
Ye 


possible to insure that y? is great at pleasure except perhaps in sub-intervals 
whose length decreases indefinitely with ry Since in the interval c+y7, d—y 


the curve is concave upward for positive y and concave downward for negative 
y, there can not be more than one zero. The discussion maybe separated 
into two parts according as y(x) has a zero or not. : 

If y(y) =0, then without loss of generality it may be assumed that within 
the interval y, d—yz, y is positive and hence, since y,,>0, that y, is positive. 
Multiplying the inequality (19) by 2py, and integrating we get 


Jf 2py.(py.) de> M’ 2ypy.de, 
yd Y 
and hence : 


p’ (x) y(")= p(x) yz (%) — p(y) yi(y) 2M’ [min p] Jf 2yy.da 
v 
=M’'y’(x) min p, (20) 





* This is at once evident if the equation is taken in the form Yoo + (9 + dK) y= to which it may be 
reduced by a transformation. 
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for all x in the interval y, d—y. Similar reasoning establishes the same 
result in the interval c+~z,y. Hence, for the case that y vanishes, we have 


2 eee 
throughout the interval c+yz, d—y the inequality i > = eo which as 





noted above, is sufficient to establish the theorem. 

If on the other hand y has no zero in the interval, let us take the function 
positive and denote by y the point at which it takes on its minimum. It may 
be assumed that y is not at the right-hand end of the interval; that special 
ease could be treated in an analogous fashion. By multiplying y by a con- 
stant, y(y) can be made equal 1, while the sign of the integral D(y, ¢, d) is 
not altered. Since y,(y)20, we have y>1 in the neighborhood of y; in fact, 
for any arbitrarily small but fixed y we can, by taking MW’ large enough, have 
y(y+n)>2. For, since y21 it follows from the inequality (19) that 


py. py.—Py.(y) =f (py2)de2M'f yde=M' (a—y). 
¥ v 
Integrating again, 
eg x mn E M’ (w—y)? 
(y—1) max p=max p f y.de> f py,dx = M f lay )dexe =. 
fh J st vy iad 


- nit >1 we have y—1>1 or y>2. From 
2 max p 

this it follows that on multiplying the inequality (19) by 2py, and integrating 
we get by a process similar to that used in (20) the formula 





By taking x—y 27 and choosing 


p(x) y2(x) >’ [ 2ypy.de > M’[y?(«)—1] min p> = y’ (x) min p. 
v 


Hence in the interval y+», d—yz we have’ - a — And since in the 
interval c+, y—yz a similar inequality may be obtained, this completes the 
proof of the lemma. 

It follows in the same way that if k>0, except perhaps at the end points 
of the interval c, d, then 4 can be taken so large negatively that for solutions 
of (15) the integral D(y, c, d) is positive. 

We have then proved by these Jemmas that in any sub-interval in which k 
has one sign, D(y) can be made positive by taking 4 large enough positively 
and also by taking 4 large enough negatively. It follows that D/(y, 0, 1) 
would be positive, and from the formula 


al al 
af kytdx = | (py2—qy’) dz=D(y, 0,1) >0 
0 0 
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and Theorem IV, that for |A| large enough the zeros of Y move to the left 
with increase of |A|. We can then enunciate 

Tuerorem VII. There exists an integer nz such that for n=n, there are 
precisely two solutions of (15) (16) oscillating n times. 

If proper restrictions are imposed on G(a#,A), a theorem similar to the 
preceding may be proved for the general equation 

(DYz)2+G(%, A)y=0. (10) 

The proof used in Theorems VI and VII, and Lemmas I and II is valid also 
under the following hypotheses on a function G’(z, 4) obtained by subtracting 
from G(z,4) a function q(x) + [G’ (x, 2) =G(a, A) —q(zx) ]. 

(a) In one or more sub-intervals of 0, 1, G’(x, A) is a monotone increasing 
function of 4 such that lim G’(#,4)20 and for at least a part of each sub- 


A= 
interval lim G(x, 2) =+o. 


A\=a 
(8) In the remaining sub-interval or sub-intervals of 0,1, G’(z,4) isa 
decreasing function of % such that lim G’(a, 2)<0, and for at least a part of 


\=—00 
each sub-interval lim G’ (2%, A) =—o. 


A\=—0 


(y) For values 4’, A” of A there are , zeros of the solution Y within the 
interval 0,1 and for a>2’, 


os = (A) G’ (a, a), and for 1<A”, 2 <—$,(A)G' (a, A), 
where ¢,, 9, are positive functions defined in the intervals 4’, +0 and —a, 2” 
respectively. 
Tueorem VIII. Under the conditions (a) (8) (y.) there exists an integer 
n, such that for n>n, there are precisely two solutions of (10) (16) oscillating 
n times. 
When the parameter value and corresponding solution of (15) (16) are 
complex, let us set A=o+it7, y=u-+iv. The differential equation resolves 
itself into the two following: 
(pu,),+tqu+sku—tkv=0, u(0)=u(1)=0, (21) 
(pv,),+quvu+ckv+tku=0,. v(0) =v(1) =0. (22 
Corresponding to A=c—it there is then a solution y=u—iv. On multiplying 
(21) and (22) in the first place by wu and v, and in the second by v and —u 
respectively, adding and integrating, one obtains the two formulae 


f i (w2 + v2) —q(u?+v*) |d~=0, fk +0%) da=0. (23) 











302 RuicHnarpson: Contributions to the Study of Oscillation Properties 


More generally one obtains by this process the formula 
B 
[p(u,v—v,u) ee | k(u?-+v®) da. (24) 


The orthogonality of two solutions ({(kY, Y,,dx=0, m+n) holds as well 
for complex as for real characteristic numbers A,,4,,. Hence by separation 
into real and imaginary parts we get the 

Tuoeorem IX. If Y,,(%,An)=UntWn, Y,=u,t+iv, are two solutions of 
(15) (16), then 

SF (qty ars) 42 =0, S Fundy +40) de=0. 


Theorems of other types may be derived of which the following is an 
example: 

TureorEmM X. If in the interval 0,1, the function k(x) changes sign once 
only, then the roots of the real and imaginary parts u,, v, of the solution 
Y,=vu,+1v, of (15) (16) separate one another. 


d (“) 
: ae (k(w+ v°) der. 
0 


da vy. 


For, formula (24) may be written 


p(uzo—v yl) = fk (u+0*) da, or p 


From the hypothesis that k changes sign but once it follows that the integral 


d (+) 
ean not vanish within the interval, and since a >0O except at the end 
dx 


f ie 
points, > sa monotone function and the theorem follows at once. We may 


note further that under the hypotheses of the theorem neither w nor v can 
gain or lose a zero at a point within the interval. For, were there a zero lost 
or gained at x=a, both y(a) and y,(a) would vanish, and formula (24) could 
be written f°k(u°’-+-v’)dx=0, which would give a contradiction. 


§5. Reduction of the General Boundary Conditions to Normal Forms. 


In the preceding sections the discussion has dealt mainly with the simple 
boundary conditions y(0)y,(0)=y(1)y,(1) =0. To facilitate the discussion 
of the most general boundary conditions it is desirable to obtain normal 
forms.* The linear self-adjoint equation of the second order will be taken in 
sii (x(x) u,(2)),+T (a, 2)u(z)=0, m(x) >0, (25) 








* The classification used in this section follows that of Hilbert and Haupt, loc. cit. The geometrical 
form into which the transformation is thrown was suggested by my colleague, Prof. H. P. Manning. 
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and if the linearly independent boundary conditions 


au (0) + au, (0) + agu (1) + au, (1) =0, 9 
2,0 (0) + Bu, (0) +6,u(1) +8,u,(1) =0, (29) 
are to be self-adjoint, the condition 
u(0) u,(0) u(1) u,(1) 
*()) F0) 3%, (0) |=") | aay (a) ai, 








must be imposed, where w%, u are any functions with continuous derivatives 
which satisfy (26). We have the relation 


| u(0) %,(0) 


Ay Ae Az U4 
Bi Bs Bs Bs 


as may be seen by applying the usual rule for multiplying determinants and 
noting from (26) that each of the four elements of one product determinant 
is the negative of the corresponding element of the other. Let us denote by 
B,; the determinant formed by taking the 7-th and j-th columns of the matrix 


(28) 


u(0) u,(0) u(1) %,(1) 














| u(1) w,(1) 





t ~ @ & a 
|G: @ Bs B 


Considering the two determinants of the w’s as the variables, we have from 
the theory of linear equations that a necessary and sufficient condition for the 
solution of (27) and (28) is that 


(1) B,—72(0) Bs,=0. (29) 








Hence B,,. and By are simultaneously zero or different from zero. 
Let us subject the dependent variable to the transformation 


u(a)=xn(a)y(x), n(x), (30) 


by which the number of zeros of the solution remains unaltered. The equation 
(25) takes on the self-adjoint form 


(PYz)et+G(%,2)y=0, pH=my'>0, G(x,a)=nyn.,+antln’, (31) 
and the conditions (26) are replaced by a similar set 


719(0) +yey.(0) +ysy(1) +ys.y,(1) =9, \ 


> 32 
8 y(0) + &rye(0) +8y(1) +8iye (1) =0. | _ 
We have at once from (31) the important formula 
1 1 
or wdae= aG ydx, (33) 


0 OA 0 OA 
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Writing u,, U2, Us, uy, for u(0), u,(0), u(1), w,(1), and using a similar 
notation for the 7’s and y’s we can interpret our problem in relation to the 
tetraedron of reference as the investigation of the line (26) subject to the 
condition (27) or (29). If B.,=0, the us and u, can be eliminated from (26) 
at the same time, giving 

Bit, + Boytz=0 or Byuwy—Byu,=0, (34) 
either of these equations representing the plane determined by the given line 
and the edge 12(u,=0, u.=0). Hence the line (26) intersects this edge of the 
tetraedron. But if B,,=0 we have from (29) B,,=0, and the straight line 
intersects also the edge 34. The equation of the plane determined by the 
given line and the edge 34 can be written 

Byglg t+ Byug=0 or Bosutg—Byut.=V. (35) 
From (34) and (35) we see that under the hypothesis B,,.—B;,=0 the boun- 
dary conditions may be written 


Case I. ou, tu,=0, tu3+u=0, (36) 
h 
—— pa Bu Bu 5 Bu __ Bos 
Bog By’ By By 


The parameters o and t may have zero or infinite values when the line (26) 
lies in a face of the tetraedron. 

In general, when any two £6’s with complementary indices are zero, the 
line intersects two opposite edges of the tetraedron and the conditions (26) 
may be reduced to a normal form similar to (36). 

The transformation (30) takes the form 

U=MYiy Us=NYitMY2, Us=NYs, Us=MYstnNsYs, m=O, "20, (37) 
and corresponds to a rotation of the face 2 of the tetraedron about the edge 
12, and of the face 4 about the edge 34, leaving the faces 1 and 2 and the 
edges 12, 13 and 34 unchanged. If we denote by d,; the determinants of the 
matrix of the coefficients y;, 6; of (32) we can read their values from tlie 
identity 
MA, +N, 1%, NzAgtMsAy Ns hy 


73 3s Ya Me il . 
18, +728, nm Po nsPstniGs nse F 








thus A,=7iBy., Ass=73By,, ete. As may be seen from these formulae and 
(31), we have corresponding to (27) or (29), the relation 
p(1)Ap—p(0)Ay=0 (38) 


which, when added to (32), makes these boundary conditions self-adjoint. 
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The discussion of the boundary conditions under the transformation may 
be sub-divided as follows: 

I. If the given line intersects the edges 12 and 34 (By=B,=0), the 
same will be true after the transformation since A,,=A,=0. To reduce the 
condition to the normal form (36) (which we shall in general use) the trans- 
formation is superfluous. This case is called the sturmian. 

From the formulae (37) it follows at once that the conditions (36) can by 
proper choice of 7;, 72, 73, 7, be reduced to one of the simpler forms 

(a) w=ys=0; (b) m=%=0; (Cc) w=ys=0; (d) y=y,=0, (39) 
these corresponding geometrically to the cases where the given line coincides 
respectively with the edge 13, the new edge 14, the new edge 23 or the new 
edge 24. 

II. If the line does not intersect the edges 12, 34, but does intersect 13, 
the transformation may be so determined that it will also intersect 24; e. g., 
by making the plane 2 pass through the intersection of the given line with the 


plane 4. In this case since u, is unaltered we have 7;=1, y,=0, and hence to 
equate A,,=73(7:5:1;+725.,) to zero we need only to choose 2 =a, After 
1 23 

the transformation we have A,,=A,,=0, and by elimination in (32) we get by 
aid of (38) | 

Case II. Yi=hys, hp(0)y,=p(1)y%, 
where 
_ As Au PO) Aw _ p(1) Arg 


h =— = — ; 
Ay, © Ay p(0) dy, p(0) Ass 





III. If the line does not intersect any of the three edges 12, 34, 13, 
(By 0, By -=0, By =—0) we can determine the transformation so that the 
face 2 of the tetraedron shall pass through the intersection of the line and the 
face 3, while the face 4 passes through the intersection of the line and the 
face 1. Thus, in the new tetraedron, the line intersects the edges 23 and 14, 
and we have the relations Ay=73(7By4+72B4) =0, As=m (usBo+.Bu) =0, 





which determine the ratios a = desired in this case. The equations of the 
line can now be written : 
Case III. w=lp(1)y, lwW(0)y=—ys, 
where 
_ Ay _ Ag ni Aj, A Ay 


~ p(l)Ay, p(l)Ais p(0)Ay = (0) Ag” 
39 
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THerorEM XI. By a change of dependent variable the self-adjoint equation 
(25) with self-adjoint boundary conditions (26) (27) may be reduced to an 
equation (31) which is again self-adjoint, and for which the self-adjoint boun- 
dary condition may be written in one of three forms characterized respectively 
by Ay=Ay=0; AxysdAy=0; Ayp=Ay=0: 

Case I. cy(0)+y,(0)=0, ty(1)+y,(1) =0, o, t constants (including o), 
Case II. y(0)=hy(1), hp(0)y,(0) =p(1)y,(1), h=constant, 
Case III. y(0) =/lp(1)y,(1), lp (0) y,(0) =—y(1), l =constant, 
The number of zeros of the solutions of the equation remains unaltered under 
this transformation. 
The value of the integral u’da is the same as that of the corre- 


ia ‘OG .. 
sponding integral { a, y dx. 
Corotuary I. In Case I the transformation may be so chosen that the 
boundary conditions can be written in one of the special forms 
(a) y(0) =y(1)=0; (b) y(0) =y,(1)=0; 
(¢) y,(0)=y(1)=0; (d) y,(0)=y,(1) =0. 


§6. Oscillation Theorems for the Sturmian Boundary Conditions. 

We propose now to study oscillation theorems for solutions of the differ- 
ential equation under the first or sturmian case of the general boundary 
conditions (26) (27). As was shown in the preceding section this case, char- 
acterized by the relations B,,=B,,=0 (A,;,,=Aj,=0), may be reduced to a 
study of the equation 


(PYz)z+ G(x, A)y=D), (41) 
under the boundary conditions 
sy(0)+y,(0)=0, ty(1)+y,(1) =0, (42) 


where co and 7 are constants, the important special forms 
(a) y (0)=y(1)=0; (b) y (0) =y,(1) =0; 


(c) y.(0)=y(1)=0; (d) y,(0)=y,(1) =0; (43) 


being considered as included for the values 0 and o. It was further shown 
that the boundary conditions may always be reduced to one of these special 


forms. 
We proceed now to a study of the nature of the dependence of 4 on co and 7. 


If 7 and y are defined as the following fundamental solutions of (41), 
y (0, A) =0, y,(0, A)=1; y(0, A) =l1, y,(9, a) =0, (44) 
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any solution y(z) can be written y=c,y+c, . For determination of c, and ¢, 
we have on substitution in (42), 


oe+¢e,=0, t[qy(1,a)+ey(1,4)]+ay,(1, 2) +e%,(1, A) =0. 
A necessary and sufficient condition that there be values of c,, c, satisfying the 


equations is that 
1 o 


tY(1,4)+Y9.(1,4) ty (1,4) +y.(1, a) 


Hence each of the parameters o, t is a one-valued function of the other and 
of a. If we assume that a zero of the solution moves continuously to the 
right or left, and note what happens as two consecutive zeros pass through 
the point ~=0, we see geometrically that for any given o the condition 
oy(0)+y,(0)=0 must be satisfied at some stage of the process. It is possible 
by analysis to ascertain precisely what happens when o or t passes over the 
range —o, +o. Corresponding to c=—o, 0, +o one has respectively 
y(0) =0, y,(0) =0, y(0) =0, and we shall show that under proper conditions 
one and only one zero of y(x) is lost or gained by the process. A similar 
change in t has a corresponding result. 
If in formula (13) we substitute from (42) and the formulae 
on da dy,(0 Oy (1 dy,(1 
0 + 00) +  =0, 2B + y(t) + Se =0, 


obtained from (42) by differentiation with regard to 2, we obtain the relation 





al Y 


0G , lo . 
(a. yda=—p(0) Gey *(0) + p(1) & = y'(1). (45) 


2’O 


dz 
Hither of the parameters o, t may be held fixed. If 7 is fixed - 7A =0, and in 


order that 4 and o be one-valued functions of one another over py oo —o, 
+oo of o it is then necessary and sufficient that the integral (57 = y'dau have 
one sign. From such considerations one may deduce a variety E: ar. of 


which the following are examples. 


TueorEM XIII. Jf 2},, 4),., denote two successive characteristic numbers 
of a solution of (41) for the boundary conditions 43(a), then in order that 
there be for all o an intermediate value of 2% corresponding to which there is 
precisely one solution of (41) for the boundary conditions cy(0)+y,(0)=0, 


* OG 
y(1)=0, it is necessary and sufficient that the integral J Oa y’dx have one 


sign for all solutions y(a) concerned, 
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THroreM XIV. If An, Anyi denote two successive characteristic numbers 
of a solution of (41) for the boundary conditions 
oy(0)+y,(0)=0, y(1)=0, c=constant, 
then in order that there be for all < an intermediate value of 4 corresponding 
to which there is precisely one solution of (41) (42), it 1s necessary and 


1 
sufficient that the integral f OF yaw have one sign for all solutions y(x) 


concerned, 
THeoremM XV. If corresponding to a',A” there are solutions of (41) (48a) 
1 
with n, and nz zeros respectively, and if f oF de has one sign, then in the 
0 
interval 4’, 2” there is one and only one value of A for which (41) (42) has a 
solution y(x) with n zeros (n<n<M,). 
We have seen in $5 that by a transformation of the dependent variable 
the equation (41) may be thrown into the form 


ats il ~ OG _, “OG , * 

(py.).+ Gy=0, cA oa Y dx =f Oa Y dz, p(x) >0, 
where the boundary conditions (42) assume one of the special forms 
Y (a1) =Y (42) =0; Y(o1) =F, (42) =0; YY, (41) =Y (42) =03 Y, (01) =Y, (a2) =0. 
As an extension of Theorems IV and IVA we have then the following: 

TureorEM XVI. If a zero of one of the functions cy+y,, tyt+y, ts held 

fixed, then with increasing 4 the zeros of the other move closer to the fixed 
zero or further away according asf 2 yan is positive or negative, the inte- 


gration extending over the interval between the zeros. 


From the preceding theorems we can deduce still further results. 


wl 
THeoreEM XVII. If J SG ytd has one sign, then between any two con- 
0 


secutive 2’s corresponding to solutions of (41) (48a) or of (41) (43d) there 
will be one each corresponding to solutions of (41) (48b) and (41) (48c), and 
between two consecutive a’s for (41) (48b) or (41) (48c) there will be one 
each for (41) (48a) and (41) (48d). 

_Corotuary: Denoting by 2, and A, two adjacent characteristic values for 
solutions of the problem (41) (48a), (41) (48d) respectively, and by 2,, A, 
and 24, a, the next greater and next smaller sets for solutions of the problem 
(41) (43b), (41) (48c) respectively, then if the number of zeros of the solution 
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corresponding to A, 1s m, the number corresponding to A, is m—1. When the 
integral is positive the number of zeros corresponding to A, and A, is m, and 
to A, and A, is m—1, while if it 1s negative, the reverse is true. 


If we consider the special form* of (41) where = > 0 and Hypothesis A 


of § 2 is satisfied while G(a, —#«) <0 and =—o in at least some portion of 
the interval, we can trace the various values of 4 corresponding to successive 
suites of ¢ and 7, and state 


TurorEM XVIII. There exists under these hypotheses an infinite set of 


characteristic numbers 2,<A,<A3...., with limit point at +o only, for 
which exist solutions of (41) (42), the solution y,(x) corresponding to 
A, (n=1, 2, 3,....) having n—1 zeros within the interval. 


The results for the orthogonal case ($1, IL) of the special equation 


(pyz)2+ (a+ak)y=0, (46) 
are contained in the foregoing theorem. In the other cases it is preferable to 
investigate the zeros through a discussion of those of the equation obtained by 
the transformation used in proving Theorem XVI. We have 


(DY.)2+(G+ak)¥=0, p(x) >0, ¥(0)¥,(0) =¥(1)¥,(1) =0, 
al al 1 
py.—qy’ — ky C= 2 a 
J, (py:— dy ) da J ydx J ky dx 


If the new equation is of the polar form we see at once that there are 
precisely two solutions of (46) (42) with n zeros in the interval. If the equa- 
tion is of the non-definite form there are two integers ny, n., such that for 
n<n, there are no solutions of (46) (42), for n>n, there are at least two, 
while for n>n, there are precisely two (Theorems VII; VII A). 


§7. Oscillation Theorems for Case II of the Boundary Conditions. 


We shall next consider the exceptional case By,.=a,3,—a,3,—0 for wh ich, 
as has been seen in § 5, the boundary conditions may be written 


y(0)=hy(1), hp(0)y,(0)=p(1)y.(1), h=constant$0. (47) 
Defining the particular solutions 7(x, 2), y(z, A) by (44), and substituting the 
solution y=c,y+c.y in (47), we obtain equations for ¢, and c, 
e=h[ay(1,a)+ey(1,a)], hp(0)q=p(1) [ay.(1, 4) +ey,(1, 4) ]. (48) 
A necessary condition for a solution is the equation 
hy(1,a hy(1,a)—1 
= y(1, a) y(1, a) ae (49) 
p(1)y,(1,4)—hp(0) p(1)y,(1, a) 





* This is a slightly more general condition than that imposed by Birkhoff. Loc. cit. 
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Since for the two solutions we have the well-known formula 
p(x) [y,(&@, A)Y (x, A) —Y, (x, A)y (x, A) ] =constant=—p(0), (50) 
this condition reduces to 
h®p (0) y (1,4) —2hp (0) + p(1)y,(1, 4) =9, 


from which we get by again using (50) 


y(1,2) Vp(0)y(A, a) 
So long as 7,(1,4) and ¥(1,4) have opposite signs there will be two values of 
h for each value of 2; when they have the same sign there will be none. We 
may distinguish three critical cases for which the solutions pass from real to 





h 





complex: 

(a) If (1,2) =0, y,(1, 4) #0, the equation (49) can by means of (50) 
be written [hy (1, 2) —1]*p(1)¥,(1, 2) =0, and since y and y¥, can not vanish 
together, it follows that hy(1,~)—1=0. Hence every element of the deter- 
minant D vanishes except that in the lower right-hand corner; it follows then 
that c.=0. The solution y(#)=c,y(1,4) vanishes at x=0, and since h#0, 
we have from (47) that it vanishes also at y=1. The function y(z) is there- 
fore a solution of (41) (43a). 

(@) If (1,4) 0, ¥,(1, 2) =0, we have in the same way c,=0 and the 
solution y=c,y also a solution of (41) (43d). 

(y) If y(1, a) =0, y,(1, 4) =0 it follows immediately by the same reason- 
ing that.all elements of D vanish. Hence c, and c, may take on any values 
and 4 is a double characteristic number. 

On multiplying equation (41) by y and integrating under the boundary 
conditions (47) we obtain the relation 


vl 1 
J py2dx = Gy'dx. (52) 
0 0 


By substitution in (13) from formulae (47) and the further formulae, 





dy(1) __ dy(0) = dh dh | dy,(0) _ dy.(1) 9. 
hae 5 Ys Gy PO) ye (0) +hp (0) “3, =p(1) “aa, (53) 


obtained from (47) by differentiation with regard to 4, we obtain the funda- 
mental formula 
' OG 1 0G 1 OG 
] eee» P h? ~ pda ater 
dh YS, ma oo ey a J, a (54) 





da 2p (0) y(0) (0) 2p(1)y(0)y.(1) 200) yA) y(0) 
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Let us denote by 4, a characteristic number of our problem for h=0; this 
is a solution of the sturmian case (41) (43b). Near it, as we know from 
Theorem XVII and its corollary, there is another characteristic number A, 
corresponding to h=o and to (41) (43c), the two solutions for 2, and 2, 
having the same number of zeros. The parameter value 2, may be equal to, 
greater than, or less than ”,. Of the aggregate of characteristic numbers for 
the two problems (41) (48a), (41) (48d) let us denote by 4,, the greatest of 
those smaller than 2, (or 4,) and by 44, the smallest of those greater. Within 
the interval 2,,, 4,2 it follows from Theorem XVII that there is but one zero 
and one infinity of h; in other words, but one solution for each of the prob- 


al 
lems (41) (43b), (41) (48¢). On the hypothesis that { oG y'dx has one sign 
0 


dA 
we are now in a position to prove that in this interval of 2 the function h(A) 
defined by (51) is monotone on both of its branches (the one including h=0 
and the other h=o). 

To prove this let us in the first place note that for h=0 or h=o we have 


a sturmian case for which existence theorems have already been established, 


and then show that as h passes through either of these values, . does not 
change sign. We see at once from (47) that y,(1,4(0)) =0; hence y(1,4(0))=0, 


and it follows from (47) that y(0,A(h)) changes sign with h. <A reference to 


the first part of (54) shows that = will then retain its sign. In the same 


way y(1,4(«0))=0, y,(1,4(0))#0; y,(0,A(h)) changes sign as h goes 


ee th be ; : ; ; 
through infinity and a retains its sign. This establishes the result since it 


al 
follows from (54) that while { SO yd retains one sign the only possibility 
0 


eee | 
of =: changing sign is when h or y(0) or y,(0) changes sign. 


By definition y(0,2(0)) =0, y,(1, 2(0)) =0, and without loss of generality 
it may be assumed that y,(0,2(0))>0. The discussion then divides itself 
into two parts according to 


Hyporuesis I. y(1,4(0))>0. Hypornesis Il. y(1,4(0)) <0. 


Under the first hypothesis the number of zeros of y (including that at 
x=0) is even, while under the other it is odd. Roughly speaking, we shall see 
that a gain or loss of a zero comes when h goes through 0 or o. Concerning 


the function 2(h) there is now sufficient data to sketch the graph. We shall 
vl 


first discuss the problem for the assumption { = ydx>0 and Hypothesis I. 


0 
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It has been proved above that = can not change sign on either branch of 


the function. On the branch containing h=0, A=2, we see from (54) that 
dh 


7A, <0, since for that particular point y(1)>0, y,(0)>0. For h=w,a~=A,, 


the number of zeros is the same as that for h=0 (Corollary, Theorem XVII) 
and hence under Hypothesis I, y,(1) and y(0) have opposite signs. From the 


second part of (54) we see then that = >0 on this other branch. Since at the 


ends of the interval h is not 0, we know that * in cases (a) and (@) either y(0) 


or y,(0) will vanish according as 4,,, 24, belong to the problem (41) (43a) or 
(41) (43d); it follows then from (54) that = oe 

Under Hypothesis II it is readily shown by the same processes that in an 
interval which we shall call 4,,, 2, to distinguish it from the other, the branch 
of the function 2 (h), which contains h=0, is monotone increasing and the branch 


containing =o is monotone decreasing. As 4 increases this form of curve will 


1 
always alternate with that obtained under Hypothesis I. When f oe y’dz <0, 
0 


the second form of curve occurs under Hypothesis I, and the first under 
Hypothesis II. 

Since by definition the intervals can not overlap, the curves can not 
overlap. If the upper bound of one interval is the lower bound of the next, 
we have a double point for A [cf. (y) above], and one branch of one curve 
unites with one of the other to form a function monotone throughout. This is 
what takes place, for example, in the case of the solutions of the equation + 


Yr2 tay=0. 





*In case (y), y(0) and y,(0) may be chosen arbitrarily; cf. next succeeding foot-note. 


1 a = 
} For this special equation y= Vix sin Via, y =cos Va, and equation (49) becomes 


2h —(h?+1) cos Vi=0, 


from which we obtain the formula h=sec VA+tan VX. For every positive \ the function h is double- 
valued. The interval \aa,X’ca iS (2m)*x*, (2m + 1)2x*, while at (2m + 1/2)°x°, h becomes zero on one branch 
and infinite on the other. The monotone decreasing branch h,=sec Vx—tan Vx may be considered as 
joined at both ends to monotone decreasing branches in the next intervals. These intervals (2m—1)°z’, 
(2m) 2x? and (2m + 1)2x*, (2m + 2)*n* are of the type Aaa A’aa- The function h, decreases from +0 to—o 
in the interval (2m —1/2)*n’*, (2m + 3/2)°x*. The function h.==sec Vx + tan Vx increases monotonely 








from —o to +o in the interval (m+ 1/2)*n°, (m+ 5/2)°x*, the branches of adjacent intervals uniting 
as in the other case. Each of the curves h,; and h, cuts two of the other set orthogonally, the points of 
intersection occurring at the end points m*n° of the interval daa, N'aae 
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It is now easy to write down oscillation theorems for this case. We note 
in the first place that under Hypothesis I the number of zeros is increased by 
unity as h goes through zero, and under Hypothesis II it is decreased by unity 
so that for h>0 the number of zeros is always even, and for h<0 it is odd. 
There are then two solutions with an even number of zeros in the one case and 
two with an odd number in the other. 


THEOREM XIX. If in an interval of % there exist two integers m,, mz 
positive or zero such that there are solutions of 


(PYz)2+G(x, A)y=0 (55) 


and (48a) with m, and m, zeros respectively within the interval 0,1, then 
1 
under the hypothesis that { SG yds has one sign throughout, there are for 
. 0 
h>0 two solutions of (55) (47) when m is even, and none when m is odd 
(mmm): for h<0 there are two solutions when m is odd aud none when 
m is even. . 
: 0G } 

The special case where Oa, >0 and G becomes negative for A=—», and 
moreover negatively infinite in at least a part of the interval and G becomes 
positively infinite for A=-+o in at least a part of the interval, includes the 
case discussed by Birkhoff and the detailed theorems derived by him hoid also 
here.* 


The orthogonal problem ($1, II) of the equation 
(PY.)2+ (d+Aak) y=0 (96) 
is contained in the special case just discussed. In discussing the polar case 
we note that it follows from the special case of the formula (52), 


1 1 
0<  (pyi—ay')da=2 { ky*da, (57) 
0 0 


1 al : 
that = ydxz= { ky’dx has the same sign as A. Hence A=0 is not included 
0 “0 


in the range of values (unless one includes the solution y=0). 

THEOREM XX. In the polar case of equation (56) there are precisely two 
solutions satisfying the boundary conditions (47) and oscillating n times 
(sweG, F, S,... of 

In the non-definite case we have from Theorem VII that for sufficiently 
large values of |%| the integral on the left of (57) is positive and hence we 
can state 











* Case II corresponds to III, p. 269 in Birkhofi’s article, loc. cit. 
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Tueorem XXI. For the non-definite case of equation (56) there exist 
two integers ny, N, (N22, 20) such that for n<n, there is no solution of 
(56) (47) with n zeros in the interval; for n>n, there are at least two, while 
for n= n, there are precisely two. 


§8. Oscillation Theorems for Case III of the Boundary Conditions. 
There remains the most general of the three normal forms obtained for 
the boundary conditions, viz.: 
y(0)=Ip(1)y-(1), Ip (0)y.(0)=—y(1), T=constant0. (58) 
Since the discussion follows the same lines as that of the preceding section it 
will be abbreviated. If we define the two particular solutions 7(2, 2), 7(a, A) 
as in (44), we get in place of formulae (49), (51), 
, |e), 4) wd)y,(1,4)—1 
Fao _ V9.(1, a)Y(A, a) 
p(1)y.(1,4) ~~ Vp(0) p(1)%,(1, 2) J 














The critical values are when 7,(1, 2) =0 and y(1, 2) =0 and the cases may be 
classified as before. 
(a) If (1,4) =0, 7,(1,4)#0, equation (59) can be written 
(1, a) (lp (1) 9, (1,4)—1]2=0 or Ip(1)¥,(1, a) =1. 


It follows that c,=0 and the solution y(x,2)=c.y(%,4) is a solution of the 
equation with the sturmian boundary condition (43¢). 

(8) If ¥(1,2)#0, ¥,(1,4)=0 we have c,=0 and y=c,y, which is a 
solution of (41) (48b). 

(y) If 9(1,2)=0, y,(1, 2) =0, it follows as before that D’ vanishes 
identically and 4 is a double parameter value. 

To replace formula (53) and (54) we have 


dy (0 il oy,(1 
oy) = & p(t)y.(1) +p (1) Ye, 


OA di 
dy,(0) . dl dy (1 
ip(0) HO) 4 F n(oyy,(0) =— HY, 
“9G ,, »f'aG., “a@ , 
dl 1, a var PY, a ede he 





eee, ae 





dz 2p(0)y,(0)y(0)  2y(1)y(0) ~ 2p(0)y,(0)p(1)y, (1) 
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It is readily shown as in $7 that = does not change sign as / goes through 


the values 0 or 0. Let us denote by A,, 2, two adjacent characteristic num- 
bers for the cases ]=0 and =o, respectively (in other words for (41) (43a), 
(41) (48d)), by A,, the greatest of the aggregate of characteristic numbers for 
(41) (48b) and (41) (48c) which are smaller than 4,, and by 4,, the smallest 
of the aggregate larger. It is readily shown as in the previous section that 
within the interval ,,, 2;,., / is a monotone function on the branch containing 
1=0 and on the branch containing h=o. Let us consider first the case 1=0; 
then y(0,2(0)) =0, y(1,2(0)) =0 and we can assume y,(0,4(0))>0. There 
will be two cases to distinguish according as we make 


Hypotuesis I. y,(1,2(0))>0; Hyporuesis II. y,(1,4(0)) <0. 


Under Hypothesis I the number of zeros is always odd, and when / goes 
through zero from negative to positive it may be’seen from (58) that two 
zeros of the solution y(z) are lost. Under Hypothesis II the number of zeros 
is even, and, as / increases through zero, two zeros of y(x) are gained. It fol- 


1 
lows from the last part of (60) that the integral { oe y’dx and = have 


dn 
opposite signs. This fixes the sign of s on one branch. To fix the sign of 


i on the other branch let us consider 4(«0)=A,. Under Hypothesis I we can 


argue from the corollary to Theorem XVII that y(0, 2(0)) and y(1, A(w)) have 

; ‘ 1 0G dl (co) 

th sign, and f the third part of (60) that wh aaa Basie. ees J 
e same Sign, and from ird part of (60) that w en f a9, Y de>, rr 

is positive. Under Hypothesis II we can easily prove that the situation is 

dl 

> da 


‘ 
other. When { be y’dx <0 the two types of curves are interchanged. With 


reversed being positive on the branch through h=0, and negative on the 


increasing 4 we see then that in both cases curves of much the same form as 
the two varieties in the preceding section alternate with one another. If 
adjacent intervals have the same end-point, a branch from the one will unite 
with a branch of the other.* 

From these data various theorems may be deduced, of which the following 


is typical. 





* For the special equation yz2-+ Ay=0 one can set up the formula for /(A) as in foot-note, p. 312. 
It is found that 1—=— ese Vi + cot Vx/ VX and that these monotone functions extend from —« to+o, 
each of the one set cutting two of the other set in double points \= (n + 1/2)?x°* of 1 (A). 
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TurorEM XXII. If in an interval of 4 there exist two imtegers u,< Ms, 
positive or zero, such that there are solutions of (55) (43a) with uw, and uw, 
zeros respectively within the interval 0, 1, then 


1 
(1) Under the hypothesis that { a y’dx be positive, there are, when l 
0 


is positive, two solutions of (55) (58) which have 2m or 2m—1 zeros in the 
interval (uy,<2m—1<2m< uu.) and there are, when | is negative, two solutions 
which have 2m or 2m-+1 zeros (u<2m<2m+1<u). 


‘ | 
(2) Under the hypothesis that { a y’dx be negative, there are, when l 


is positive, two solutions of (55) (58) which have 2m or 2m+1 zeros in the 
interval (uy.<2m<2m+1<4u,) and there are, when | is negative, two solutions 
which have 2m or 2m—1 zeros (u,)<2m—1<2m< 4). 


One can make the results of this theorem more specific by giving the con- 
ditions necessary to characterize the branch of the function /(A) which is 
involved. Birkhoff has done this for the special case treated by him (loc. cit., 
p. 269, I, IL), but we shall content ourselves with stating that the same classi- 
fication may be made in the general case treated here. 

Theorems for the special equation (56) analogous to Theorems XX, XXI, 
can be at once written down. 


Brown UNIVERSITY, January, 1918. 





